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1.0  INTRODUCTION 


1.1  Philosophy 

This  year's  research  effort  concentrated  on  evaluating  the 
state-of-the-art  in  material  modeling  for  dynamic  loading  conditions  and 
developing  some  basic  tools  for  further  research. 

Because  complex  dynamic  loading  of  soils  occurs  under  many  different 
situations,  it  has  been  necessary  to  gather  information  from  various 
disciplines  within  soil  mechanics.  The  areas  of  earthquake  and  ground-shock 
seem  to  dominate,  but  other  areas  such  as  wave  and  wind  effects,  vehicle 
loading,  machine  vibrations  and  impact  loading  also  are  very  much  concerned 
with  dynamic  soil  properties.  The  magnitude  and  nature  of  forces  imposed  on 
tne  soil  varies  quite  widely  among  these  types  of  loading.  As  a  result,  the 
people  working  in  a  specific  area  tend  to  concentrate  on  the  soil  behaviors 
wnich  seem  to  be  most  important  for  their  problem.  The  basic  philosophy 
behind  this  research,  however,  is  that  a  soil  mass  really  has  only  a  single 
fundamental  behavior.  The  properties  of  a  soil  mass  cannot  be  different 
just  because  the  load  is  a  vibrating  machine  or  a  100  kt  nuclear  blast. 
Rather,  certain  properties  will  dominate  soil  behavior  under  certain  load 
conditions.  A  truly  accurate  and  complete  constitutive  model  would  be  based 
on  the  most  fundamental  properties  and  cover  all  behaviors.  Unfortunately, 
in  ..^eloping  a  comprehensive  soil  model,  one  is  limited  by  testing 
capabilities,  both  in  the  lab  and  field,  by  the  ability  to  mathematically 
express  and  use  complex  models,  and,  to  a  varying  degree,  one  is  limited  by 
imagination  in  recognizing  fundamental  soil  properties. 

In  light  of  this,  our  effort  is  directed  toward  gaining  a  better 
understanding  of  both  soil  constitutive  models  and  insitu  soil  behavior.  In 
engineering,  one  must  obtain  the  right  mixture  of  accuracy  and  simplicity. 
If  a  soil  model  which  is  fundamentally  correct  can  be  suitably 
compartmentalized,  and  if  one  knows  what  properties  are  most  pertinent  to 
the  problem  at  hand,  then  the  necessary  elements  of  the  model  can  be  chosen 
and  the  soil  behavior  may  be  suitably  modeled. 
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1.2  Scope 


This  report  deals  primarily  with  material  modeling  in  its  current  state 
as  it  applies  to  soils.  A  general  discussion  of  soil  models  and  their  use 
is  included.  A  description  of  a  computer  code  which  is  presently  being 
developed  under  this  effort  and  which  will  be  used  to  study  material  models 
is  presented.  In  addition,  there  is  a  section  on  insitu  material  properties 
summarizing  a  series  of  two-dimensional  finite  difference  calculations  which 
has  been  done.  And  finally,  some  special  topics  which  are  important  to 
material  modeling,  shear  behavior  and  effective  stress  modeling,  are 
di scussed. 


1.3  Methodology 

In  order  to  proceed  with  minimum  duplication  of  effort,  a  review  of  the 
existing  literature  was  undertaken.  The  general  topics  covered  were: 

i)  Material  Modeling 

a)  Review  Papers 

b)  Specific  Soil  Models 

c)  Calculational  Studies/Modeling 

ii)  Insitu  Testing 

a)  Explosive  Loading 

b)  Insitu  Shear  Properties 

iii)  Laboratory  Testing 

Appendix  A  is  a  summary  list  (by  topic)  of  the  more  important 
references.  Please  note  that  this  list  is  continuously  being  updated  and  is 
certainly  not  exhaustive. 

Evaluation  of  various  material  models  has  begun,  using  a  computer  code 
which  is  capable  of  simulating  several  different  stress  and  strain  loading 
paths.  The  computer  program  allows  model  study  without  having  to  perform 
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expensive  boundary  value  calculations.  It  is  by  no  means,  however,  a 
replacement  for  finite  diff erence/ element  studies  because  it  cannot 
reproduce  overall  geometry  and  material  interaction  effects.  In  view  of 
this,  several  finite  difference  calculations  of  Cylindrical  In-Situ  Test 
(CIST)  sites  were  also  performed.  This  allowed  the  study  of  material 
behavior  under  insitu,  large  strain  conditions.  Because  of  the  sizable  CIST 
data  base,  the  calculations  also  served  as  a  basis  for  evaluating  a  modeling 
effort  with  actual  field  data. 

This  report  summarizes  the  past  years  efforts  as  outlined  in  the 
statement  of  work  contained  in  the  following  section.  Soil  modeling  is 
reviewed  in  Section  2,  the  Soil  Element  Model  is  described  in  Section  3  and 
in-situ  material  behavior  is  presented  in  Section  4.  Section  5  contains  tne 
results  of  two  special  study  areas  and  the  effort  is  summarized  in  Section  6. 

1.4  Statement  of  work 

The  following  are  the  tasks  as  outlined  in  Contract  No. 
F49620-80-C-0088  "Fundamental  Properties  of  Soils  for  Complex  Dynamic 
Loadings". 

TASK  1.  Development  of  Plotting  Routines  for  Data  Display 

ARA  will  prepare  plot  routines  for  2-D  finite  difference  codes  to 
allow  plots  of  the  states  of  stress  and  strain,  particle  velocities 
and  displacements.  Included  will  be  the  variation  of  these 
parameters  at  various  times  as  well  as  plots  of  the  CIST  data. 

TASK  2.  Perform  Two-Dimensional  Finite  Difference  Calculations  of  the  Middle 
Gust  Met  Site  CIST  (CIST-9)  and  Dry  Lake  Valley  CIST  (CIST-23) 

Using  the  best  available  cavity  pressure  description  for  CIST-9  and 
CIST-23,  ARA  will  perform  the  2-D  calculations  using  the  following 
material  models: 

•  Elastic  Based  on  Seismic  Wave  Speeds 
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•  Engineering  Model  (AFWL)  Fit  to  Laboratory  Recommended 
Properties 

•  Engineering  Mooel  Developed  from  1-D  Analysis  of  CIST-9  by  AFWL 

•  Hydrostatic 

•  Elastic,  Perfectly  Plastic  Using  Seismic  Wave  Speeds  and 
Laboratory  Failure  Surface 

•  Engineering  Model  with  Revised  Properties  to  Better  Fit 
Experimental  Velocity-Time  Histories 

•  Engineering  Model  Updated  (on  basis  of  preceding  calculation) 
to  Better  Fit  Data 

TASK  3.  Perform  Two-Dimensional  Finite  Difference  Calculations  of  CIST 
Events  in  Saturated  Sand 


A  series  of  calculations  as  described  for  CIST-9  will  also  be 
performed  for  saturated  sands.  CIST-10  and  CIST-15,  16  will  receive 
the  primary  focus.  Perform  additional  calculations  as  necessary  to 
uetter  resolve  the  oasic  response  of  the  dry  clay,  saturated  clay, 
clay  shale,  dry  sand,  and  saturated  sand  which  will  be  studied  in 
Tasks  2  and  3.  It  is  estimated  that  8-10  additional  calculations 
will  be  required. 

TASK  4.  Interpretation  of  Calculational  Studies  to  Describe  the  Fundamental 
Response  of  the  Two  Geologic  Materials 

The  calculations  performed  in  Tasks  2  and  3  will  be  used  to  develop 
an  understanding  of  tne  effect  of  the  assumption  concerning  the  form 
of  the  material  model  on  the  behavior  of  tne  materials.  The 
difference  exhibited  in  the  stress  and  strain  paths  by  the  various 
models  will  be  compared  to  the  implied  behavior  of  the  real 
materials.  Comparisons  of  these  stress  paths  to  those  utilized  in 
laboratory  tests  and  expected  to  occur  in  field  conditions  will  be 
used  to  evaluate  the  potential  of  the  models  studied  and  to  suggest 
improved  laboratory  testing  procedures.  In  particular,  the 

characteristics  of  constitutive  models  that  would  be  required  to 
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reproduce  the  measured  response  of  these  two  materials  will  be 
described  and  compared  to  other  available  models  that  were  not 
evaluated  in  this  study. 


Recommendations  concerning  constitutive  model  development  necessary 
to  better  model  these  materials  will  be  developed. 

TASK  5.  Development  of  a  Soil  Element  Code 


A  computer  code  will  be  developed  to  allow  study  of  the 
characteristics  of  the  various  material  models  to  prescribed  stress 
or  strain  paths.  The  code  will  be  general  enough  to  reproduce  the 
stress  paths  representative  of  current  laboratory  tests  as  well  as 
paths  of  interest  which  cannot  be  reproduced  in  current  test 
apparatus. 

This  code  will  then  be  exercised  for  representative  material  models 
to  develop  a  better  understanding  of  the  limitations  and  advantages 
of  current  models  and  laboratory  testing  devices.  At  least  one 
model  from  each  of  the  following  groups  will  be  evaluated: 

•  Elastic 

•  Viscoelastic 

•  One-Dimensional  Curve  Fit 

•  Ideal  Elastic-Plastic 

•  Elastic  Non-Ideally  Plastic 

As  a  minimum  each  model  evaluated  will  be  subjected  to  the  following 
loading  conditions: 

•  Repetitive  Pure  Shear 

•  Hydrostatic  Compression 

•  Uniaxial  Strain 

•  8iaxial  Compression  (Std.  Triaxial  Test) 

•  Plane  Strain  with  Rotation  of  Principal  Strain  Direction 
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•  True  Triaxial  Loading 

•  Biaxial  Extension 

•  Complex  Loading  Histories 

In  all  cases  the  response  will  be  evaluated  for  various  levels  of 
loading  and  repetitions  of  load.  The  complex  loading  histories  will 
be  developed  from  study  of  CIST,  earthquake,  and  explosive  surface 
burst  calculations. 

TASK  6.  Evaluation  of  Shear  Response  in  In-Situ  Material  Property  Tests 
ARA  will: 

Perform  theoretical  and  calculational  studies  of  the  shear  stresses 
generated  in  the  plane  (DISK  HEST),  cylindrical  (CIST)  and  spherical 
dynamic,  high  stress,  in-situ  tests.  These  studies  will  focus  on 
the  nature  (location  and  magnitude)  of  the  shear  stresses  and 
strains  developed,  their  spacial  and  temporal  variation,  and  how 
they  might  be  measured  and  interpreted. 

Analyze  experimental  data  from  representative  tests  of  each  type  to 
corroborate  the  theoretical  studies  and  establish  experimental 
measurement  limitations  on  the  recovery  of  shear  response  data. 

Develop  alternatives  to  current  in-situ  test  configurations  and 
measurements  to  improve  recovery  of  shear  data. 

TASK  7.  Theoretical  Analysis  of  Special  Problems  in  Modeling  Soil  Behavior 

Study  the  consequences,  in  terms  of  computer  storage  requirements, 
of  using  true  multi-phase  models  in  finite  difference  codes 
especially  with  respect  to  the  AFWL  CRAY  computer. 

Evaluate  alternate  approximate  methods  of  treating  effective  stress 
in  finite  difference  calculations. 
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Study  the  physical  basis  for  soils  increase  in  strength  at  very  high 
(several  kilobars)  pressures. 

Other  issues  identified  in  previous  tasks. 

The  general  result  of  this  task  will  be  input  to  future  theoretical 
work  on  model  development. 

► 

» 

> 

* 

» 

> 

» 

I 
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2.0  SOIL  MODELING 


2.1  Theoretical  Requirements 

Constitutive  relations  used  to  model  physical  phenomena  must  obey  the 
laws  of  thermodynamics  and  other  principles  of  continuum  mechanics.  (A 
model  may  consider  the  soil  material  to  be  particulate,  but  homogeneity  is 
usually  still  assumed).  The  model  should  produce  behavior  which  is  smooth 
and  continuous,  as  this  seems  to  be  the  case  for  soil,  and  infinitesmal 
changes  in  load  should  produce  infinitesimal  changes  in  state. 

Several  material  models  in  use  today  do  not  satisfy  all  the  above 
criteria.  This  is  primarily  a  result  of  the  "ad  hoc"  nature  of  soil 
modeling  to  date.  There  is  a  tendency  to  develop  a  model  which  produces  the 
required  behavior  under  the  given  conditions  with  little  concern  for 
generality,  completeness,  or  theoretical  correctness.  These  models  are  then 
usually  applied  to  different  conditions  by  an  unsuspecting  modeler  who  soon 
discovers  some  unusual  behavior. 

2.2  Behavioral  Requirements 

Certain  behavioral  characteristics  of  soils  have  been  observed  in  lab 
tests,  insitu  tests,  and  (in  practice)  under  actual  field  loading 
conditions.  Although  it  is  impossible  to  completely  summarize  soil  behavior 
in  this  report,  a  complete  soil  model  would  need  to  address  each  of  the 
following: 

a)  Soils  exhibit  non-linear  response,  even  at  strains  as  low  as 

1(T5. 

d)  Permanent  (irrecoverable)  deformations  occur  under  both  pressure 
and  deviatoric  loading. 

c)  Initial  unloading  (shear  and  volumetric)  is  very  often  elastic  in 
nature. 
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There  usually  is  a  convex  (or  "stiffening")  shape  to  the 
volumetric  pressure-strain  response  while  there  is  a  concave  (or 
"softening")  shape  of  the  shear  stress-strain  behavior. 

Soils  display  a  shear  strength,  t  v,  beyond  which  shear 

max 

"failure"  takes  place  -  usually  defined  as  the  development  of 
large  deformation  under  small  increases  in  applied  stress. 

Normal  stress  (confining  pressure)  has  a  significant  effect  on 
properties. 

There  is  a  coupling  of  volume  behavior  and  shear  behavior  often 
referred  to  as  shear-di latancy  coupling. 

The  response  of  a  soil  up  to  a  particular  stress  state  depends  on 
the  stress  path  taken  to  oDtain  that  state  (stress  path 
dependency) . 

Soils  may  exhibit  anisotropy  due  to  their  history,  i.e. 
deposition,  overconsolidation,  compaction,  etc. 

Stress  induced  anisotropy  can  occur  as  a  result  of  current 
loading  conditions. 

At  low  stress  levels,  strain  will  occur  in  the  direction  of  a 
stress  increment,  while  at  high  stress  levels,  strain  occurs  in 
the  direction  of  the  overall  stress  state. 

Soils  may  display  a  "peak"  strength  with  subsequent  strain 
softening  to  a  "residual"  strength. 

Cyclic  loads  produce  special  effects  such  as  "ratcheting"  of 
hysteresis  loops.  If  the  cycled  stress  level  is  low  enough  an 
eventual  equilibrium  state  will  be  reached.  Energy  dissipation 
(hysteresis)  occurs  even  at  low  strain  levels. 

Multi-directional  cycling  produces  different  effects  than 

unidirectional  cycling. 
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o)  The  properties  of  soils  vary  with  strain  magnitude. 

p)  Soils  may  exhibit  loss  of  structure  (sensitivity). 

q)  Deformation  in  some  soils  may  become  highly  localized,  e. g. 
formation  of  "shear  bands"  in  clay  soils  under  triaxial  test 
conditions. 

r)  Rate  of  deformation  affects  some  soil  properties. 

If  one  considers  that  soil  is  actually  a  three-phase  material 
(air-water-soil  solids)  there  are  many  additional  behaviors  which 
must  be  accounted  for.  Some  of  the  more  important: 

s)  Liquefaction  due  to  excessive  pore  pressure  and  subsequent  loss 
of  shear  strength. 

t)  Pore  fluids  migrate  depending  on  a  number  of  gradient  fields. 

u)  Pore  fluids  create  additional  rate  effects. 

As  of  yet,  no  soil  model  has  taken  all  of  these  behaviors  into  account, 
and  this  may,  in  fact,  not  be  possible. 

2.3  Practical  Requirements 

To  successfully  apply  a  particular  constitutive  model  to  a  practical 
engineering  situation,  it  must  meet  certain  requirements.  First,  it  must  be 
understandable.  A  model  which  is  so  complex  to  the  point  of  being  overly 
cumbersome  is  not  very  useful  to  an  engineer.  Second,  the  model  must  be 
physically  meaningful.  If  each  element  of  the  model  can  be  associated  with 
one  or  more  of  the  above  soil  behaviors,  the  model  is  much  easier  to  work 
with  and  alter  as  necessary.  Finally,  a  model  must  be  suitable  for  use  on 
existing  computer  hardware.  This  is  a  somewhat  limiting  factor,  but 
capabilities  are  improving  with  each  new  generation  of  computer,  and  this 
greatly  broadens  the  alternatives  for  choosing  a  suitable  model. 
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2.4  Available  Soil  Models 

The  following  is  a  list  of  the  major  models  which  are  availaole  for 
soil.  (Each  of  these  may  have  a  number  of  subsets.)  The  models  which  have 
been  installed  in  the  soil  element  are  summarized  in  Section  3.3. 
Information  on  the  others  may  be  found  by  consulting  the  appropriate 
references  as  outlined  in  Appendix  A. 

Simple  Models 

a)  Elastic  (Linear  or  Non-Linear) 

b)  Visco-Elastic 

c)  Elastic-Perfectly  Plastic 

d)  Curve-Fit  Models 

More  Complex  Models 

e)  AFWL  Engineering  Model 

f)  Cam-Clay  Model 

g)  Cap  Model 

h)  Rate  Type  Fluid  Models 

i)  Lade's  Model 

j)  Prevost's  Model 

k)  Endochronic  Models 

l)  Multi-phase  Models 
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3.0  SOIL  ELEMENT  MODEL 


3.1  Purpose 

There  are  many  constitutive  models  for  soils  which  have  been  developed 
and  tested  to  varying  degrees.  Some  models  have  their  basis  in  elastic 
theory  or  the  theory  of  plasticity.  Others  are  developed  specifically  to 
produce  one  or  more  of  the  behaviors  discussed  earlier  (Section  2.2, 
Behavioral  Requirements).  In  order  to  bring  several  models  into  direct 
comparison,  a  computer  code  has  been  written  which  can  drive  these  models  at 
the  user's  discretion.  This  "Soil  Element  Model",  is  a  tool  for  performing 
material  modeling  studies.  It  has  been  developed  with  generality  in  mind  so 
that  it  may  be  useful  in  several  ways.  Parameter  studies  of  a  particular 
model  may  be  performed,  which  provide  a  valuable  understanding  of  how  a 
model's  behavior  depends  on  the  assigned  coefficients.  One  may  subject  two 
or  more  models  to  identical  boundary  conditions  and  compare  the  results. 
Model  behavior  can  also  be  compared  with  test  data,  either  in  an  effort  to 
fit  a  model  to  a  set  of  data  or  to  predict  soil  response  for  a  particular 
test.  And  most  important,  the  Soil  Element  Model  (SEM)  is  useful  for  model 
development  and  syntnesis;  new  relationships  may  be  tested  or  model  elements 
may  be  combined  to  produce  the  best  model  behavior. 

3.2  Description 

The  SEM  is  a  small  computer  program  which  has  been  developed  at  ARA  and 
presently  is  operational  on  the  CRAY-1S  computer  at  the  Air  Force  Weapons 
Laboratory.  It  is  basically  an  adaptation  and  generalization  of  a  code 
called  "SANDR"  which  was  written  by  G.  Baladi  at  WES  for  Sandia  Laboratories 
(Ref.  1).  SANDR  is  capable  of  simulating  the  behavior  of  the 
effective-stress  cap  model  under  standard  triaxial  test  boundary 
conditions.  SEM  incorporates  several  material  model  options  and  can 
simulate,  at  present,  true-triaxlal,  pure  shear,  and  arbitrary  boundary 
conditions  as  well  as  standard  triaxial. 

The  program  consists  of  a  driver,  a  choice  of  constitutive  model 
algorithms,  and  plotting  routines.  The  basic  idea  is  illustrated  in  Figure 
3.1.  The  Driver's  role  is  to  handle  input  and  output  operations  and  provide 
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Figure  3.1.  Basic  soil  element  model  logic. 


2 


tne  material  model  with  a  strain  increment,  while  enforcing  any  specific 
boundary  conditions  wnicn  may  be  appropriate.  For  example,  in  the  standard 
triaxial  test,  the  strain  increment  tensor  must  be  controlled  so  as  to 
produce  a  stress  state  with  =  o-j-.  The  material  models  provide  a 
constitutive  relationship  ([H])  which  generates  a  new  stress  tensor: 


[do]  =  [H]  [de] 


(3.1) 


It  is  important  to  keep  the  material  model  distinct  from  the  ariving 
mechanism  to  allow  easy  alteration  and/or  substitution. 

All  of  the  models  are  "strain  controlled",  as  in  equation  3.1  above. 
This  is  consistent  with  most  finite  element  and  finite  difference  code 
formulations.  As  a  result,  making  a  model  follow  a  particular  strain  path 
is  relatively  easy.  Most  laboratory  tests,  however,  are  run  under  stress 
controlled  conditions.  In  studying  a  particular  model  it  would  be 

convenient  to  subject  it  to  a  specific  load  path,  i.e.  that  of  a  particular 
laboratory  test.  This  is  done  in  SEM  via  iteration.  That  is,  if  the  stress 
state  which  is  calculated  with  a  particular  material  model  ooes  not  meet  the 
boundary  stress  conditions,  the  strain  increment  tensor  is  adjusted  and  the 
stress  is  recalculated.  This  involves  finding  an  efficient  method  for 
adjusting  the  strain  increment  tensor.  The  more  general  test  boundary 
conditions,  such  as  true  triaxial,  are  somewhat  more  difficiult  to 

iteratively  duplicate. 

SEM  uses  any  set  of  consistent  units.  Output  appears  in  the  same  units 
as  are  input.  A  number  of  plot  options  may  be  chosen  to  graphically  display 
the  results. 

3.3  Constitutive  Relations 

The  soil  element  model  has  great  potential  for  aiding  in  model  study 
and  is  constantly  developing.  The  following  paragraphs  will  describe  the 
present  status  of  the  SEM. 


» 
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The  following  five  constitutive  models  are  presently  incorporated  into 
the  the  program: 

1)  Elastic 

2)  Elastic  -  Perfectly  Plastic 

3)  One-Dimensional  Curve  Fit 

4)  AFWL  Engineering  Model 

5)  Effective  Stress  Cap  Model 

Appendix  B  describes  the  notation  used  in  the  model  subroutines  and 
gives  some  general  definitions.  An  attempt  was  made  to  be  consistent  in 
going  from  one  model  to  the  next,  although  this  is  not  always  possible. 
Note  that  all  models  assume  the  material  to  be  homogeneous  and  isotropic. 

3.3.1.  Elastic  Model 

The  constitutive  relationship  here  is  simply  Hooke's  Law.  Only  two 
constants  are  required  to  completely  define  a  linearly  elastic,  isotropic 
material,  we  use  K  and  G,  where 


K  = 


G 


(3.2) 


and  S.j  ■  and 
expression  for 


e^  are  deviatoric  stress  and  strain,  respectively, 
stress  in  terms  of  strain  may  then  be  derived: 


An 


S 


ij 


ij 


°kk  U  ) 


C1j 


6kk  (s  ) 

~r  ( 'r 
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combining , 


and  substituting, 

°ij  =  2  G  eij  +  K  €kk  ^5ij^ 

•ij  -  «  S  I'lj  -  ^  +  K  lkk  (4ij> 

•ij  '  2  G  ‘M  +  [-■  K  ~  ?  °]  'kk  (‘Ij*  <3'3> 

Equation  (3.3)  describes  both  the  volumetric  and  deviatortic  response 
of  an  elastic  material. 

3.3.2.  Elastic  -  Perfectly  Plastic 

In  the  elastic-plastic  case,  a  failure  surface  is  defined  (in  Vv- 
space)  up  to  which  a  material  behaves  linearly  elastically.  When  the 
stress  path  intersects  the  failure  envelope,  plastic  strain  occurs  according 
to  an  associated  flow  rule  and  the  devatoric  stress  may  not  exceed  the 
failure  surface. 

The  failure  surface  is  defined  as  an  exponential  (see  also  figure  3.2a), 


f  (Jj) 


(3.4) 


Because  it  is  defined  in  terms  of  Vv  ,  the  intersection  of  the 
failure  surface  with  the  deviatoric  plane  is  a  circle  (figure  3.2b). 

When  a  stress  increment  based  on  eleastic  behavior  violates  the  failure 
surface  (see  figure  3.2c)  some  plastic  strain  occurs  and  a  correction  must 
be  made  to  bring  the  stress  state  back  to  the  failure  surface.  An 
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a) 


Failure  Surface 


c)  Correction  Procedure  for  Violation  of  Failure  Surface 


Figure  3.2.  Elastic-plastic  model . 


6 


associative  flow  rule  is  enforced  and  corrections  are  made  normal  to  the 
failure  surface.  First  a  volumetric  plastic  strain  is  calculated  as  follows: 


(3.5) 


(3.6) 


A  ratio  is  calculated  based  on  the  amount  of  necessary  correction,  and 
the  stress  deviators  are  multiplied  oy  this  adjustment: 


S 


ij 


f(^) 


(3.7) 


3.3.3.  One-Dimensional  Curve  Fit 

Curve-fitting  is  an  empirical  method  for  modeling  the  behavior  of  a 
material.  A  trend  is  observed  in  some  laboratory  test  data  and  a 
mathematical  curve  is  fit  to  it  which  reproduces  that  data.  The  curve  may 
be  exponential,  hyperbolic,  polynomial,  or  any  other  convenient  functional 
form. 


The  curve-fit  routine  in  the  SEM  is  aimed  at  reproducing  the  cyclic 
shearing  behavior  observed  in  soils.  It  is  taken  from  a  paper  entitled 
"Nonlinear  Soil  Models  for  Irregular  Cyclic  Loadings"  by  Robert  Pyke  (Ref. 
2).  The  soil  model  fits  a  hyperbolic  expression  to  the  intitial  loading 
shear  stress-strain  curve,  then  uses  extended  "Masing"  rules  to  model  the 
hystresis  loops.  These  rules  postulate  that  the  shear  modulus  on  each 
loading  reversal  assumes  a  value  equal  to  the  initial  tangent  modulus  for 
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the  initial  loading  curve  and  that  the  shape  of  the  unloading  or  reloading 
curves  is  the  same  as  the  initial  curve,  except  that  the  scale  is  enlargrd 
by  some  factor. 

The  model  is  basically  a  Davidenkov  -  class  model,  giving  the  shear 
stress  in  terms  of  shear  strain: 


t  = 


+  6 

max 


(3.8) 


where  n  =  1  for  initial  loading. 


and  n 


+  =  reloading 
-  =  unloading 


The  terms  of  the  above  equation  are  defined  in  figure  3.3. 

The  volumetric  behavior  is  treated  as  linearly  elastic  with 

°kk  =  3  K  Ykk 

although  this  can  be  modified  so  that  K  is  any  functional  form,  e.g. 

K  .  K  <„kk). 

3.3.4.  AFWL  Engineering  Model 

A  non-energy-dependent  version  of  the  AFWL  engineering  model  has  been 
installed  into  the  SEM.  The  model  is  elastic-ideally  plastic  with  volume 
hysteresis.  The  model  is  described  in  terms  of  a  failure  surface  (in 
Vv  vs  P  space)  and  a  hydrostat  (in  P  vs  v  space). 

Pressure  -  volume  behavior  is  controlled  by  the  hydrostat,  which  is  a 
series  of  linear  segments  (see  figure  3.4b).  Upon  unloading  the  model 
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Figure  3.3.  Hyperbolic  shear  model. 


von  Mises  Limit 


a)  Failure  Surface 


Deviatoric  Plastic  Potential 


i 

I  ..  .  .  .  —  -  _ — 

c)  Plastic  Potentials 

Figure  3.4.  AFWL  engineering  model. 
20 


allows  irrecoverable  volume  strain  as  defined  by  the  difference  in  the 
loading  and  unloading  moduli i.  Note  that  the  excess  compression,  u,  is  used 
instead  of  ekk,  the  volume  strain,  to  define  the  volume  pressure 
relation.  As  a  result,  it  is  not  quite  correct  to  call  the  hydrostat  slopes 
bulk  modulii.  Specifically: 


K  = 


wnere 


But 


M 


Thus 


—  <ekk  +  1) 
“  ekk 


€kk 


and,  BKL  =  K  +  P 


(3.9) 


(3.10) 


(3.11) 


(3.12) 


where  BKL  is  the  hydrostat  slope  and  K  is  bulk  modulus.  The  relationship 
between  BKL  and  K  is  pressure  dependent  with  the  difference  becoming 
significant  only  at  higher  pressures. 

The  failure  surface  (see  figure  3.4a)  consists  of  a  Drucker  -  Prager 
frictional  segment  in  combination  with  a  von  Mises  cut-off.  Below  the 
failure  surface,  the  material  behaves  elastically: 

s/  - 3  G  euE  * 3  *  *kkE  <V  <3-13> 

where 

K  -  K  (w) 

G  »  3/2  (K)  Q  ~ 
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At  the  failure  surface,  there  are  two  plastic  potentials  which 
effectively  uncouple  distortional  and  volumetric  plastic  strain  (see  figure 
3.4c).  The  volumetric  flow  rule  is  associative: 


d 


(3.14) 


and  the  deviatoric  flow  rule  is  non-associative,  with  stress  being 
vertically  corrected  back  to  the  failure  surface: 


d  V 


(3.15) 


d 

where,  \  =  V  (3.16) 

26  Vv 


3.3.5.  Effective  Stress  Cap  Model 

The  SEM  contains  a  version  of  the  cap  model  which  is  intact  from 
Baladi's  original  SANDR  code.  At  present  only  one  modification  has  been 
made  to  allow  for  a  curved  failure  surface. 

The  cap  model  is  an  elastic-plastic  constitutive  model  which  has  been 
developed  in  accordance  with  the  rules  of  classical  plasticity  (see,  for 
example.  Hill  (Ref.  3)).  As  a  result,  its  mathematical  description  is  a 
little  more  complicated  than  that  of  the  AFWL  model,  which  was  developed 
with  calculational  results  in  mind.  The  cap  model  was  an  important 
development,  because  it  effectively  couples  shear  and  volumetric  behavior. 

The  behavior  of  the  model  in  the  elastic  range  is  governed  by  the 
elastic  bulk  and  shear  modulii.  The  bulk  modulus  is  a  function  of  only: 
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where  K^,  and  are  material  constants  and  Gr  is  the  overburden 
stress  due  to  gravity.  G,  the  shear  modulus,  is  taken  to  be  a  function  of 
the  second  invariant  of  the  stress  deviator  tensor, 


G  -  _ [1  -  Gj  exp  {-  G2  ^/jp)]  (3.18) 

1  -  G1 

where  G^,  Gj  and  G2  are  material  constants.  The  parameters  for  the 
bulk  and  shear  modulii  are  evaluated  by  fitting  curves  to  match  the 
unloading  hydrostat  Dehavior  and  the  unloading  stress  difference  -  strain 
difference  triaxial  behavior,  respectively. 

The  plastic  behavior  of  the  model  is  described  in  stress  space  by  two 
functions  (see  figure  3.5a).  A  modified  Drucker-Prager  failure  surface  is 
defined  by: 


f^j)  -  (C$  +  C)  -  C  exp  (BJj)  +  «  (3.19) 

where  Cs,  C,  B  and  o  are  material  parameters  which  can  be  oDtained  from 
triaxial  test  data.  Cs  represents  a  cohesive  strength,  a  a  frictional 
strength,  8  defines  a  gradual  increase  in  strength,  and  C  a  maximum  shear 
strength.  By  properly  defining  the  constants,  a  curved  or  linear  failure 
surface  may  be  specified. 

An  elliptical  strain-hardening  surface  is  defined  by: 


F  (0,,  *)  -  1  {[*(<)  -  L(<c)3?  -  [J,  -  L(k)]2)  0’5  (3.20) 

JF 
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where  X(<)  and  L(tc)  are  as  shown  in  figure  3.5a.  R  is  the  ratio  of  the 
major  to  minor  axes  of  the  elliptical  yield  surface.  is  the  hardening 
parameter  and  controls  cap  movement: 


with  D  and  W  material  constants.  Thus,  hardening  of  the  cap  is  a  function 
of  plastic  volumetric  strain. 

The  cap  model  is  an  "effective”  stress  model  in  that  it  has  both 
drained  and  undrained  responses  for  the  volumetric  behavior  (figure  3.5b). 
The  undrained  model,  obtained  from  an  undrained  test,  simulates  total  stress 
behavior.  The  Grained  model,  ODtained  from  a  drained  set  of  tests, 
simulates  effective  stress  behavior.  Since  total  stress  equals  effective 
stress  plus  pore  pressure,  the  pore  pressure  response  of  the  material  is 
quite  readily  calculated. 

The  computation  of  stress  given  a  strain  increment  in  the  cap  model 
subroutine  procedes  by  first  assuming  an  elastic  trial.  Checks  are  then 
made  to  see  if  any  control  surfaces  have  been  exceeded.  If  the  failure 
surface  is  exceeded  plastic  strain  is  calculated  according  to  an  associative 
flow  rule  and  the  hardening  parameter  is  updated  to  account  for  plastic 
volumetric  strain.  Iteration  is  needed  to  maintain  consistency  between 
location  of  the  stress  point,  the  cap,  and  the  value  of  <,  the  hardening 
parameter. 


3.4  Test  Paths 

The  SEM  model  is  capable  of  exercising  any  of  the  installed 
constitutive  relations  over  the  following  test  paths: 

i)  isotropic  compression/consolidation  (undrained/  drained) 

ii)  uniaxial  strain 
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iii)  standard  triaxial  ( o ^  =  o 

iv)  cubical  triaxial  (strain  controlled,  i.e.  rigid  platens) 

v)  simple  shear  (strain  controlled) 

vi)  arbitrary  strain  paths 

Several  of  these  paths  are  shown  in  figure  3.6  for  the  triaxial  plane. 
Note  that  one  can  also  simulate  the  results  of  proportional  loading  tests, 
i.e.  those  with  arbitrary  directions  in  the  triaxial  plane. 

3.5  Sample  SEM  Exercises 

Appendix  C  illustrates  the  use  of  the  Soil  Element  Model  with  six 
examples.  The  first  five  examples  deal  with  one  material,  McCormick  Ranch 
Sand  (Ref.  4),  for  which  several  models  have  been  preliminarily  fit.  The 
models  were  then  subjected  to  various  loading  conditions,  with  the  results 
shown  by  the  figures  in  the  appendix.  The  last  example  is  of  the  effective 
stress  cap  model,  which  has  been  fit  to  model  a  deep-sea  sediment  (Ref.  1). 
Table  3.1  summarizes  the  paramters  which  are  necessary  to  describe  each  of 
the  models. 


Further  exercises  of  this  type  will  be  done  to  study  material  model 
responses.  The  models  will  be  compared  with  data  from  both  laboratory  and 
insitu  experiments. 

3.6  Continuing  Work 

The  SEM  will  play  an  important  part  in  the  dynamic  soil  properties 
effort  which  is  continuing  at  ARA.  The  following  enhancements  and 
modifications  constitute  the  direction  of  development  for  the  computer  code 
itself.  Other  additional  areas  of  improvement  may  later  become  evident. 

i)  Model  refinement.  Enhance  the  response  of  the  models  already 
installed  in  the  code.  Examples  are  allowing  for  non-linear 
elasticity,  isotropic  hardening  or  softening  of  the  elastic 
plastic  model,  and  combining  Pyke's  shear  model  with  the 
volumetric  portion  of  various  other  models. 
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IC  -  isotropic  consolidation  ozz  =  oRR  =  o-j-j 


CTC  -  conventional  triaxial  compression  crRR 
CTE  -  conventional  triaxial  extension  o ^  = 


=  Oy y  =  constant 
constant 


°RR  ”  aTT 


TC  -  triaxial  compression  a 

(constant  mean  stress  aRR  =  aTT,  - 
test,  pure  shear 

°RR 

TE  -  triaxial  extension  oRR  =  Oyy,  - 


RR 


+ 


+  o 


ZZ 


constant 


=  constant 


UX  -  uniaxial  strain  (a  varies  with  Poisson's  ratio) 


Figure  3.6.  Strain  paths  for  some  laboratory  test  configurations 
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(Stress)-1  Failure  Surface  Defines  Amount  of  failure 

surface  curvature 


TABLE  3.1.  (Concluded) 


ii)  Additional  Models.  Install  other  constitutive  models  such  as 

Lade's  model,  Duncan  and  Chang's  hyperbolic  model,  Prevost's 
model,  endochronic  model. 

iii)  Time  Dependence.  Introduce  rate  effects  into  the  SEM  by 

associating  a  time  increment  with  each  strain  increment.  Having 
accomplished  this,  such  models  as  viscoelasticity, 
viscoplasticity,  and  rate-type  fluid  may  be  studied. 

iv)  Stress  Controlled  Driving.  This  is  important  for  simulating 

laboratory  tests  (which  are  usually  stress  controlled)  and  for 
reproducing  arbitrary  stress  paths.  Rather  than  mathematically 
inverting  the  models,  it  is  anticipated  that  stress  control  will 
be  accomplished  by  iterating  until  the  desired  stress  state  is 

satisfied  for  each  increment. 

v)  Anisotropy.  Modify  the  installed  models  to  allow  for  anisotropic 

material  response.  Also  allow  for  initial  anisotropic  (KQ) 

consolidation  prior  to  tri axial  shearing. 

vi)  Plot  Enhancement.  Increase  the  available  catalog  of  plots  as 

needed.  Additional  plots  already  planned  are  principle  stress 

space  plots,  invariant  plots,  time  domain  plots,  and  model 
comparison  (overlay)  plots. 
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4.0  1NSITU  MATERIAL  BEHAVIOR 


4.1  Introduction 

Any  attempt  to  model  insitu  soil  behavior  based  on  laboratory  test 
results  encounters  several  difficulties.  First,  there  are  the  problems 
involved  with  obtaining  representative  samples  and  minimizing  sample 

disturbance.  Secondly,  the  loading  conditions  of  conventional  lab  tests 
limit  the  type  and  number  of  stress/strain  paths  to  a  relatively  simple 
few.  As  a  result,  studying  the  response  of  soil  under  insitu  conditions 
(either  controlled  test  conditions  or  actual  field  conditions)  is  a 

necessity.  This  section  deals  with  one  type  of  dynamic  insitu  testing. 

4.2  The  Cylindrical  In-Situ  Test 

The  particular  test  wnicn  will  be  emphasized  here  is  known  as  the 
Cylindrical  In-Situ  Test  (CIST)  (Ref.  5)  and  was  developed  by  The  Air  Force 
Weapons  laDoratory  in  the  early  70' s.  A  CIST  involves  the  detonation  of  a 
PETN  explosive  in  a  vertical  cylindrical  borehole  and  the  subsequent 
measurement  of  ground  motions  in  the  soil  or  rock  around  the  cavity  (fig. 
4.1).  CIST  provides  information  on  material  response  in  a  large  strain 

regime  and  under  rapid  loading  conditions.  A  large  data  base  for  many 
geologies  and  materials  has  been  developed  since  the  first  CIST  in  1971 
(table  4.1) . 

Ideally,  a  CIST  event  produces  axisymmetric  motion  about  the  centerline 
of  the  explosive  cavity.  In  addition,  it  is  assumed  that  simultaneous 
loading  occurs  along  the  entire  cavity.  In  reality,  conditions  may  be 

somewhat  different,  depending  on  the  particular  site  conditions  and  care 
with  which  the  test  is  performed.  In  an  actual  geology  there  may  occur 
lens-like  deposits,  non-horizontal  layers,  material  non-homogenet ies,  etc, 
all  undetected  by  any  subsurface  investigation.  Any  of  these  will  tend  to 
produce  non-ax isymmetric  motions,  and  this  is  occasionally  observed  where 
there  have  been  redundant  radial  measurements. 

An  attractive  feature  of  CIST  geometry  is  that  several  soil  layers 
having  different  properties  can  be  simultaneously  tested  in  what  is 
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Explosive  Cavity 
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Figure  4.1.  CIST  cavity  with  geologic  stratification 
(From  Ref.  6,  p.  10) 


3' 


TABLE  4.1.  CIST  EVENT  SUMMARY 
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MINUTEMAN  Site  M-28,  18  Oct  74  Wet  Clay  REST  TEST  III 

North  Dakota 


CIST  Mo.  Location  Date  Geology  Associated  Ground 

Motion  Event(s) 
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theoretically  a  "uniaxial"-l ike  stress  state  at  early  times.  If  this 
assumption  is  true,  models  may  be  formulated  by  using  one-dimensional  wave 
propagation  theory  and  iteratively  matching  calculational  and  experimental 
velocity  waveforms.  Usually,  however,  the  motions  are  at  least 
two-dimensional  in  nature  especially  at  late  times  (and,  unfortunately, 
sometimes  three-dimensional) .  Two-dimensional  motion  results  from  waves 
eininating  from  the  ends  of  the  cavity  and  the  ground  surface  (fig.  4.2), 
cratering  effects,  and  from  the  effects  of  layering.  For  these  reasons,  a 
two-dimensional  axisymmetric  simulation  tends  to  produce  results  more 
closely  approximating  the  experimental  results. 

4.3  CIST  Calculations  -  Description 

There  were  several  goals  in  performing  these  calculations: 

a)  To  Detter  define  the  ground  motion  experienced  during  a  CIST 
explosion. 

o)  Study  the  effects  of  varying  material  model  parameters  on  the 
calculated  ground  motion. 

c)  Arrive  at  insitu  .  models  by  matching  calculational  and 
experimental  results. 


Two-dimensional  CIST  simulation  is  basically  a  dynamic  axisymmetric 
boundary  value  problem.  Constitutive  relationships  are  chosen  and  material 
interfaces  are  specified.  The  calculation  may  be  driven  by  several  means. 
A  pressure  boundary  may  be  specified  along  the  cylindrical  cavity  which  best 
represents  measured  experimental  pressure  time  histories.  Unfortunately, 
few  reliable  sets  of  cavity  pressure  measurements  have  been  made  for  a 
CIST.  (This  is  primarily  because  the  harsh  cavity  environment  has  usually 
destroyed  one  or  more  components  of  the  measurement  system.)  Measured 
horizontal  velocity  time  histories  at  the  3  ft.  range  gages  have  been  used 
to  drive  some  one-dimensional  calculations  (e.g..  Ref.  6).  This  method 
avoids  the  uncertainty  surrounding  cavity  pressures.  The  method  used  in 
this  study  was  to  model  the  gases  in  the  cavity  using  a  simple  y-law  gas 
relation: 
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KEY: 


Prime  indicate  relief  waves. 


-  loadinq  Wave 

-  —  Relief  Wave 

P  Compression  Waves 
S  Shear  Waves 
SP  von  Schmidt-Like  Waves 
R  Rayliegh  Waves 


2.  Major  wavefronts  from  a  cylindrical  explosion 
in  an  elastic  medium. 
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p  »  (y  -  1)  P  e 
where  P  =  Pressure 
p  =  Density 
e  =  Energy 
y  =  Coefficient 

An  initial  energy,  eQ  is  introduced  at  the  beginning  of  the  calculation 
sufficient  to  provide  the  desired  peak  pressure.  From  experience  (Ref.  7) 
the  best  representation  of  CIST  waveforms  has  been: 

Pmax  -  40'° 

pmixture  =  0.00225  9/cc 

y  =  1.25 

The  decaying  pressure  in  the  cavity  produces  calculated  ground  motion  which 
seems  to  be  qualitatively  appropriate. 

AFTON,  a  two-dimensional  finite  difference  code  (Ref.  8),  was  used  to 
perform  the  calculations.  Actually,  a  specialized  version,  known  as 
CISTAFTON,  which  was  developed  at  the  Air  Force  Weapons  Laboratory  circa 
1972,  was  used.  AFTON,  in  addition  to  being  axsymmetric  as  required  by  CIST 
geometry,  has  the  feature  of  allowing  arbitrary  (generalized)  grid  motion 
(as  opposed  to  strictly  Lagrangian  or  Eulerian).  During  a  CIST  event,  a 
crater  develops  directly  adjacent  to  the  cavity  at  the  ground  surface.  If 
grid  motion  is  restricted  to  purely  Lagrangian,  the  grid  stretches  out  in  a 
narrow  zone  in  the  cratering  region  and  causes  the  necessity  for  a  very 
small  calculational  timestep.  If  done  properly,  arbitrary  grid  motion  can 
circumvent  this  problem  . 

The  AFWL  Engineering  Model  was  used  in  these  calculations.  The  AFWL 
model  was  developed  for  this  kind  of  study  in  an  effort  to  make  parameter 
variation  and  selection  somewhat  simplified.  (This  is  why  the  hydrostat  and 
failure  surface  have  been  linearized.)  The  model  and  its  parameters,  as 
they  appear  in  CISTAFTON,  are  shown  and  defined,  in  fig.  4.3  and  Table  4.2, 
respectively.  (This  version  of  the  model  varies  somewhat  from  the  version 
which  appears  in  the  SEM.)  Note  that  the  model  contains  an  energy 
dependence  which  is  not  usually  exercised  during  a  CIST  calculation. 
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Hydrostatic 
Pressure,  P 


Figure  4.3.  CISTAFTON  AFWL  Engineering  Material  Model. 
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Variable 

Symbol 

Description 

Units 

RHOZ 

=  p0 

=  initial  material  density 

g/cc 

POISL 

=  VL 

=  loading  Poisson's  ratio 

— 

POISU 

=  vy 

=  unloading  Poisson's  ratio 

— 

CL1 

=  CLi 

=  constrained  loading  velocity 

fps 

CL2 

=  CL2 

=  constrained  loading  velocity 

fps 

CU 

=  cu 

=  seismic  velocity,  unloading 

fps 

CZ 

=  cz 

=  constrained  loading  velocity, 
corresponding  to  initial  slope  of  high 
pressure/density  portion  of  hydrostat 

fps 

CV 

-  cv 

=  Constrained  Velocity,  corresponding 
to  tail  of  unloading  hydrostat 

fps 

PI 

=  Pi 

=  hydrostatic  pressure  at  break  between 
initial  unloading  and  Cj_i  portion 

psi 

P2 

=  P2 

=  hydrostatic  pressure  at  break  between 

Cu  and  Cl2 

psi 

P3 

=  P3 

=  hydrostatic  pressure  at  initiation  of 
hign  pressure/density  portion 

psi 

AM3 

=  v3 

=  excess  compression  at  P3 

— 

ST1,  2 

=T1»T2 

=  tensile  strengths  of  peak  and  residual 
failure  surfaces,  respectively 

psi 

Yl,  Y2 

=Yl.Y2 

=  '\/^2'  ~  axis  intercepts  of  failure 
surfaces 

psi 

SI,  S2 

=Si,S2 

=  slopes  of  Drucker-Prager  portion  of 
failure  surfaces 

— 

VM1,  VM2 

=  peak  and  residual  von  Mises  limits 

psi 

AMS 

=  us 

=  material  parameter,  high  energy/ 
density  (HED) 

— 

BKM 

=  km 

=  material  parameter  (HED) 

psi 

TILA 

=  A 

=  materia]  parameter  (HED) 

— 

TILB 

=  B 

=  material  parameter  (HED) 

— 

TILE 

=  eo 

=  initial  internal  specific  energy  ergs/gm  x  10*2 

ES 

-  es 

=  material  parameter,  energy  ergs/gm  x  10*2 

TABLE  4.2.  GUIDE  TO  CISTAFTON  AFWL  MATERIAL  MODEL  INPUT  PARAMETERS 
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In  an  effort  to  evaluate  the  material  sensitivities  to  the  various 
model  parameters,  a  few  simplifications  of  the  model  were  used  for  some  of 
the  calculations.  These  included  purely  elastic  (with  no  shear  failure)  and 

hydrodynamic  (elastic,  with  no  shear  strength).  In  addition,  various 
combinations  of  elastic-plastic  parameters  were  used. 

Thirteen  calculations  have  been  performed  to  date.  Emphasis  was  placed 
on  CIST  9  and  10,  which  are  wet  clay  and  sand  sites  respectively,  and  CIST 
23,  a  dry  sand  site.  Two  calculations  were  done  for  "CIST  00"  a 

hypothetical  CIST  site.  This  generic  site  is  used  for  studying  model  and 

geometry  effects.  Table  4.3  summarizes  the  calculations  and  should  be  used 

as  a  guide  for  the  following  discussions. 


4.4  CIST  Calculations  -  Results 
4.4.1.  Cavity  Pressure 

Figures  4.4  and  4.5  summarize  the  cavity  pressure  time  histories  which 
were  generated  for  several  of  the  AFTON  calculations.  It  is  seen  that 
because  the  gases  vent  from  the  top  of  the  cavity,  the  pressures  decay  much 
more  rapidly  in  this  region  than  in  the  middle  or  bottom  regions  of  the 
cavity.  Note  that  the  pressure  near  the  top  decays  to  zero  within  about  6 
ms. 


At  depth,  decay  of  cavity  pressure  also  depends  significantly  on  the 
amount  of  cavity  expansion.  This  is  controlled  oy  the  amount  of  material 
deformation  directly  around  the  cavity.  Both  the  hydrostat  and  failure 
surface  influence  deformations.  By  comparing  the  8.5  m  and  12.2  m  depths  of 
CIST  calculations  00.1  and  00.2  (fig.  4.5)  one  can  easily  see  that  a  stiff 
hydrostat  causes  the  pressure  to  "hang  up"  much  longer.  The  effect  of  shear 
failure  can  be  seen  by  comparing  cavity  pressures  for  calculations  903  and 
904.  The  hydrodynamic  case  (904)  of  no  shear  strength  allows  more 
deformation  of  the  cavity  (and  therefore  lower  pressures)  than  the  material 
with  significant  shear  strength  (903). 

In  general,  a  two  term  exponential  of  the  form: 
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TABLE  4.3.  SUMMARY  OF  CISTAFTON  CALCULATIONS 


Calculated 
out  to  ... 

(ms) 

9*6 

7.9 

o 

• 

CO 

o 

00 

8.2 

iH 

co 

13.9 

CD 

_c 

CL 

4-> 

>> 

<D 

CD 

<V  1 — 

C 

s- 

3 

s: 

s 

0) 

s 

s: 

s: 

3  <U 

r— 

> 

i- 

r—  U 

•r— 

4-> 

•*—  fO 

cO 

3 

2 

CO 

2 

2 

2 

rO 

4- 

U.  U 

Q_ 

a. 

o 

Q_ 

a. 

a_ 

3 

O 

1 

l 

S- 

1 

1 

j 

CO 

2 

Q 

Q 

cd 

Q 

Q 

o 

hsi 

4-> 

co 

S_ 

+-> 

fO 

S- 

U 

s- 

t/>  CD 

U 

a) 

u 

O 

cO 

cO 

CO 

O  CL 

■f— 

c 

•r* 

<D 

(V 

Q> 

S-  >, 

4-> 

«r— 

4-> 

4-> 

C 

c 

tz 

t/1 

i — 

to 

CO 

•r* 

«r— 

>> 

ro 

•f- 

<o 

CO 

r— 

. - 

r— 

U 

r— 

i — * 

•r* 

L&J 

1 — 

UJ 

UJ 

GO 

CO 

CO 

-a  • 

CO  4- 

to 

JD 

c 

r—  S_ 

•r-» 

cO 

c 

c 

o 

O0 

3 

E 

—1 

o 

o 

•r- 

r— 

3:  to 

CO 

•r— 

+-> 

a> 

c 

-a 

4-> 

4-> 

a. 

u 

TJ 

o  a; 

>> 

c 

<o 

CO 

•r— 

•r— 

o 

*r**  s- 

■a 

cO 

u 

u 

U 

4-> 

s: 

40  3 

o 

•r- 

a 

00 

00  1 — 

• 

CO 

cO 

_Q 

CO  •  r— 

"O 

<J 

•r— 

•r- 

a; 

r— 

cO 

r—  CO 

=■> 

“O 

•O 

a 

UJ 

— / 

LU  n- 

or 

cO 

o 

o 

o 

s: 

s: 

• 

u  • 

*— 1 

C\J 

ro 

f-H 

CM 

ro 

<—  o 

o 

o 

o 

o 

o 

o 

o 

CO  2 

CD 

CD 

CD 

CD 

o 

o 

o 

<D 

r-H 

r— < 

fH 

~o 

a>  4-> 

“O 

co  a> 

• 

c 

3  to 

4-> 

CO 

cO 

QJ 

to 

to  c 

tO 

i — 

r— 

cO 

r-  O 

>> 

cO 

r-  QJ 

CD 

<u  -r- 

>* 

cO 

-C 

<o  > 

f— 

“O  4-> 

r~ 

CO 

CO 

C 

o  »o 

r— 

O 

o  u 

1 

s  — 

«  >> 

>> 

a> 

U>  CD 

c—  a; 

3 

•r-  rtf 

X 

CO 

CO 

as  c 

*—  O 

4->  *— 

T3 

r*— 

*— 

CO 

.  -a 

u  o 

«0  r- 

s-  o 

C 

O 

<d 

-C 

+J  c 

O  -M 

•r-  CO 

CO 

cO 

to 

CO  co 

CD  to 

«-  o 

Q-  >i 

to 

4- 

• 

CO 

c 

a> 

•O 

4- 

o 

>) 

■o 

•  *f— 

-M  L. 

»  c 

4-> 

fO 

cO 

— N  C 

4->  co 

cO  o 

CM  ns 

CD 

4-> 

L. 

r— 

C\J  co 

CO  U 

s:  •*- 

' — to 

3 

to 

Lu 

o 

' — ^to 

CO  o 

• 

• 

• 

• 

• 

• 

• 

fH 

CM 

*3* 

to 

r-C 

CM 

■o 

U 

c 

*5  c 

o 

to 

o  -* 

•r— 

•  *f— 

i  CD  ZD  CD  O 

r-4  O 

*  U 

O  -t-> 

CD  4->  -O 

■M  r—  CO 

2  cO 

(—  -I-  <T3 

t—  aj  a. 

U 

00  UJ  OO  S- 

co  2  O 

V-  o 

—1 

o 

4-> 

<0 

CD  O  +■»  *"“ 

U  C<+J 

Q  <u  o 

UJ 

3 

CD 

•  t*  <_> 

o 

3E 

to 

a> 

CD 

CO 

a. 


x 

a; 

c 

“O 

0) 

3 

c 

*r— 

+-> 

c 

o 

<d 


43  APPLIED  RE/EARCH  AXfOCIATEMOC. 


1 


CIST  No.  and  Material  Models  used  Calc.  Description  Hydrostat  Failure  Calculated 

Location  for  Calculation  Set  No.  Type  Surface  Type  out  to  .. 
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Cavity  Pressure  (MPa) 


(4.1) 


P(t) 


-c2(t) 


+  C3  e 


-c4(t) 


where  Cj,  C^,  C^  and  C4  are  constants,  determined  to  fit  the 
calculated  waveform,  is  adequate  for  mathematical ly  describing  these 
pressure-time  histories.  These  fits  could  then  be  used  as  approximate 
driving  pressure  boundaries  for  one-dimensional  calculations.  The  pressures 
generated  here  are  also  important  if  one  wishes  to  compare  experimental  data 
with  calculated  cavity  environments  in  order  to  check  validity. 
Unfortunately,  tnis  is  an  area  still  open  to  question  due  to  lack  of  high 
confidence  cavity  pressure  measurements. 

4.4.2.  Generic  CIST 

Tne  hypothetical  site  conditions  for  CIST  00  calculations  00.1  and  00.2 
are  shown  in  figure  4.6.  The  first  calculation  was  for  a  uniformly  elastic 
site  and  the  second  was  for  a  two-layer  elastic  site,  with  soft  material 
overlaying  stiff  material. 

Figure  4.7  illustrates  the  basic  response  differences  for  the  fast  and 
slow  materials.  The  stiff  material  transmits  the  stress  pulse  at  a  much 
higher  frequency  and  at  a  somewhat  higher  level,  but  sustains  much  less 
material  deformation. 

Velocities  for  the  two  cases  are  compared  in  figures  4.8  thru  4.11. 
The  frequency  oifference  is  again  readily  apparent.  Decay  of  velocity  and 
stress  is  compared  in  figures  4.12  and  4.13.  According  to  Reinhart,  (Ref. 
9),  the  theoretical  decrease  in  stress  with  range  for  a  cavity  of  infinite 
extent  would  be: 


°H  ‘  (WO1'2  (4.2) 

where  a  =  Cavity  Radius 
r  =  Radius 

Pmax  =  Peal<  Cavity  Pressure 
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gure  4.6.  Generic  CIST  calculational  profiles  (CIST  00). 


Depth  =  12.2  m 


V 


Figure  4.11.  Velocity  time  histories  for  2-layer  elastic  (00.2) 
calculation  at  depth  =  12.2  m. 


Peak  Horizontal  Velocity  (m/s) 


30 


Figure  4.12.  Variation  of  peak  horizontal  velocity  with  range  for  CIST  00 
calculations. 


5' 


This  is  also  plotted  in  figure  4.13.  The  differences  observed  between 
theory  and  calculation  are  accounted  to  the  finite  length  of  cavity  (two 
dimensional  effects)  in  a  CIST  experiment.  This  is  illustrated  by  the  fact 
that  the  theory  agrees  best  with  the  calculation  of  the  stresses  at  the  mid 
height  of  the  cavity. 

The  finite  two-dimensional  nature  of  the  CIST  geometry  can  be 
summarized  by  looking  at  the  stress-strain  behavior  (fig.  4.14).  Note  the 
vertical  strains  near  the  top  and  bottom  of  the  cavity  due  to  relief 
effects.  The  cylindrical  effects  are  shown  in  the  tangential  stress-strain 
response.  In  general,  the  stress  state  becomes  very  complex  after  passage 
of  the  initial  wavefronts. 

End  effects  and  layering  of  materials  with  different  wave  propagation 
properties  will  also  affect  site  response.  The  increased  vertical  motion 
near  the  interface  of  calculation  00.2  is  an  example  of  this  (fig.  4.9). 
Velocity  vector  plots  (figs.  4.16  and  4.16)  can  sometimes  illustrate  these 
effects.  In  figure  4.15,  for  example,  one  can  easily  see  the  shear  wave 
propagating  from  the  bottom  of  the  cavity  at  later  times.  In  figure  4.16, 
one  can  pick  out  waves  propagating  both  due  to  end  effects  and  layering 
effects.  The  shear  wave  coming  off  the  bottom  of  the  cavity  (see  4.2  ms  VV 
plot)  can  also  be  picked  of  the  shear  stress/strain  plots  or  the  principal 
stress/strain  angle  rotation  plot.  Anticipated  arrival  at  a  5  ft.  radius 
from  the  cavity  bottom  is  calculated  as  follows: 


Cs  =  (1070  m/s)  (0.53) 

Cs  =  567  (m/s)  (4.3) 


Figure  4.17  points  out  the  passage  of  this  shear  wave  as  seen  in  the  time 
histories. 
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H-Horizontal 

•  •  •  •  V-Vertical 
—  •  —  *  T-Tangential 
S-Shear 
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Figure  4.14.  Stress-strain  behavior  for  elastic  case 
at  range  a  1.5m. 


vector  snapshots  for  an  elastic  CIST. 


Figure  4.16.  Velocity  vector  snapshots  for  2-layer  elastic  CIST 
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Figure  4.17.  Shear  wave 


4.4.3  AFWL  Model  Variation 


A  series  of  calculations  for  CIST  9  was  performed  in  which  the  material 
model  type  was  varied  to  study  the  effect  on  resulting  ground  motion. 
Figure  4.18  gives  site  geology  and  layering  for  these  calculations. 

Four  cases  were  run.  Case  901  was  an  elastic  calculation,  with 
wavespeeds  determined  from  CIST  9  arrival  times  as  shown  in  figure  4.19. 
Case  902  used  models  based  on  laboratory  data,  as  determined  by  WES  (Ref. 
10).  Case  903  was  elastic,  as  in  901,  but  used  the  lab  failure  surfaces  as 
used  in  902.  Case  904  was  elastic  with  no  shear  strength  and  is  called  the 
"hydroaynamic"  case.  The  material  models  for  the  four  calculations  are 
shown  in  figures  4.20  and  4.21. 

The  results  of  the  model  study  are  summarized  in  figures  4.22  and 
4.23.  The  elastic  case  (901)  shows  a  peak  velocity  with  rebound  and  a 
subsequent  negative  phase  of  velocity.  At  the  other  extreme,  the 

hydrodynamic  material  (904)  flows  outward  due  to  lack  of  shear  strength. 
(Figures  4.24  and  4.25  show  the  resulting  large  deformations  for  this  case. 
Note  the  significant  upward  velocity  component  in  the  upper  soil  layers.) 
For  case  903,  the  upper  layers  do  not  have  appreciable  strength  to  behave 
significantly  different  from  the  hydrodynamic  material.  At  the  5.8  m  depth, 
however,  shear  strength  is  high  enough  to  bring  the  velocity  waveform  down 
quite  a  bit,  and  limit  continuing  outwar H  flow. 

As  is  expected,  based  on  the  hydrostats  the  lab  models  yield 
considerably  slower  propagation  velocities  for  the  soil  layers.  Also 
because  of  the  softer  nature  of  the  hydrostat,  the  particle  velocities  are 
much  higher.  The  material  does  not  have  much  shear  strength  and  tends  to 
maintain  outward  velocity  after  failure.  Notice  that  at  the  shallow  depth 
the  lao  models  are  in  close  agreement  with  the  hydrodynamic  model  but  at  the 
deeper  depth  they  are  more  like  the  elastic-plastic  models.  By  comparing 
results  of  type  with  data  for  the  events,  general  conlcusions  can  be  made 
with  regard  to  the  controlling  portion  of  the  material  model. 
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•  EGG  MGIII  data  report  -  preliminary  seismic 
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Figure  4.21.  Failure  surfaces  for  CIST  9,  calculations  902  and  903 
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Figure  4.25.  Velocity  vector  snapshots  for  hydrodynamic  CIST  9. 


4.4.4.  CIST  10  and  CIST  23 


CIST  10  was  conducted  in  the  saturated  coral  sands  at  the  Eniwetok 
Proving  Ground.  The  variation  of  materials  with  depth  was  relatively  minor 
and  the  water  table  at  the  time  of  the  experiment  was  at  the  ground 

surface.  The  ground  motion  data  also  showed  little  variation  with  depth  and 
quite  simple  wave  forms  at  all  locations.  The  particle  velocity  rose 
rapidly  to  a  peak  value  and  after  a  slight  recovery  began  to  indicate 

essentially  constant  horizontal  velocity  continuing  for  significant  times. 
Close  to  the  CIST  cavity  this  velocity  increased  slightly  with  time  and  at 
the  more  distant  ranges  was  uvistant  at  approximately  the  peak  value  reached 
during  the  initial  rise.  This  indicates  that  the  material  reached  the 
failure  surface  during  the  initial  loading  and  then  began  to  flow  outward 
much  like  a  heavy  fluid.  The  material  model  parameters  discussed  in  the 
previous  section  indicates  that  this  behavior  could  be  reproduced  with  a 
very  simple  model.  In  fact  the  hydrodynamic  model  yielded  wave  forms  quite 
similar  to  those  observed  in  CIST  10.  The  model  used  to  fit  these  data 

involved  a  bilinear  loading  hydrostat  with  a  small  amount  of  air  and  a 
von-Mises  type  failure  surface  with  a  low  value  of  shear  strength.  Because 
the  material  was  essentially  saturated,  a  Poisson's  ratio  of  0.49  was  used. 
This  allowed  the  material  to  load  up  below  the  failure  surface  but  fail 

almost  immediately  upon  unloading.  Figure  4.26  shows  representative  data 
and  the  calculated  waveforms  for  the  3.7  m  depth.  The  comparison  here  is 
good,  indicating  that  a  quite  simple  model  is  quite  adequate  for 
representing  the  behavior  of  the  saturated  sand.  Initially,  the  fit  of  the 
material  parameters  necessary  to  represent  this  data  was  obtained  almost 
completely  by  examination  and  knowledge  of  the  laboratory  testing  results. 
Only  a  couple  of  iterations  were  necessary  to  produce  the  agreement 
indicated  on  figure  4.26.  The  material  model  used  for  this  calculation  is 
shown  in  figure  4.27.  Clearly,  for  this  saturated  sand,  a  very  simple 
material  model  is  quite  adequate. 

This  behavior  is  contrasted  to  that  of  a  dry  sand  (CIST  23)  and  a  clay 
shale  (CIST  9)  in  figure  4.28.  The  most  obvious  differences  between  the  wet 
and  dry  sand  are  the  initial  arrival  time  and  the  rise  time  to  peak 
velocity.  The  early  time  comparisons  between  the  clay  shale  and  the 
saturated  sand  indicate  very  similar  performance,  however,  at  late  time  the 
essentially  strengthless  wet  sand  continues  to  flow  outward  whereas  the  clay 
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b)  Failure  surface  (0-4.6  m)  and  representative  stress  path. 


Figure  4.27  CIST  10  hydrostat  and  failure  surface  (calculation  1003). 


Horizontal  Velocity 


shale  tends  to  recover  and  exhibits  negative  velocities  during  the  times  of 
comparison.  This  figure  illustrates  the  large  differences  in  response 
exhibited  by  the  different  materials  for  the  same  basic  loading  condition. 
This  suggests  tnat  not  only  will  the  parameters  in  a  model  be  quite 
different  for  the  different  materials  but  that  different  response  mechanisms 
may  be  governing  the  behavior  of  the  material.  If  this  is  true,  then  these 
experiments  should  provide  a  test  for  the  generality  of  the  basic  model 
formulation. 

The  CIST  23  data  shown  on  figure  4.28  is  repeated  in  figure  4.29  and 
compared  to  the  calculated  behavior.  The  AFWL  Engineering  Model  used  here 
was  developed  based  on  laboratory  data,  and  then  improved  by  iterating  the 
solution  to  ODtain  a  better  match  to  the  experimental  data.  Notice  that 
neither  of  the  calculations  match  the  exact  qualitative  features  of  the 
experimental  data  at  all  ranges.  In  general,  however, the  model  developed  by 
iteration  is  a  reasonable  fit.  Comparisons  of  this  kind  indicate  that  for 
dry  sand  materials  a  somewhat  more  complex  model  may  be  required  in  order  to 
reproduce  Doth  the  qualitative  and  quantitative  features  of  experimental 
data.  A  similar  comparison  for  CIST  23  from  the  12.8  m  depth  is  shown  in 
figure  4.30.  In  this  case  both  the  laboratory  derived  model  and  the  model 
developed  by  iteration  are  in  quite  good  agreement  with  the  qualitative 
features  of  the  experimental  data  and  reproduced  the  magnitudes  quite  well. 
The  model  developed  by  iteration,  however,  seems  to  reproduce  the  detailed 
features  of  the  waveform  better  than  the  model  based  on  laboratory  data. 

4.5  Discussion  of  Results 

The  CISi  calculations  have  helped  to  better  isolate  the  behavior  of 
soils  under  insitu  conditions.  By  comparing  various  geologies  under  similar 
loading  conditions,  fundamental  differences  in  material  behavior  may  be 
studied.  In  addition,  the  basic  wave  types  produced  in  a  CIST  test  can  be 
better  quantified  by  calculations  of  this  type. 

The  comparisons  presented  In  the  previous  sections  have  illustrated 
that  the  cylindrical  insitu  test  also  provides  a  reasonable  method  to 
evaluate  various  material  models  and  to  check  their  generality  against  a 
reasonably  simple  insitu  experiment  of  a  complex  dynamic  nature.  The 
ability  of  a  material  model  to  reproduce  the  qualitative  nature  of  wave 
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Figure  4.29  CIST  23  velecity  waveform  comparisons,  depth  =  6. 
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forms  generated  in  the  CIST  experiments  provide  a  first  order  screening 
process  for  models  which  have  previously  demonstrated  a  capability  to 
reproduce  laboratory  data. 


It  is  clear  that  two-dimensional  CIST  calculational  studies  contribute 
to  a  better  understanding  of  the  effect  of  material  properties  on  observed 
material  response  and  should  be  pursued.  A  number  of  problems  with  the 
AFTON  code  can  be  circumvented  by  using  either  the  recent  update  of  AFTON 
(CRALE,  California  Research  Arbitrary  Lagrangian  -  Eulerian)  or  STEALTH  a 
newer  code  which  we  are  currently  using  on  other  calculational  programs. 

Additional  calculations  will  also  study  the  effects  of  model  variations  and 
will  further  examine  basic  material  response  differences.  In  this  way,  CIST 
remains  a  valuable  tool  in  studying  insitu  soil  behavior. 
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5.0  SPECIAL  PROBLEMS  IN  MODELING  SOIL  BEHAVIOR 


5.1  Snear  Behavior 

The  study  of  the  calculations  as  presented  in  the  previous  sections  has 
indicated  tnat,  while  the  CIST  provides  a  good  method  for  studying  the 
behavior  of  materials  and  the  effect  of  various  material  models  on  that 
uehavior,  the  technique  is  not  ideally  suited  for  identifying  propagating 
shear  waves  upon  which  to  base  estimates  of  the  material  shear  modulus.  The 
primary  shear  waves  in  a  CIST  experiment  are  generated  at  the  top  and  bottom 
of  tne  cavity  and  propagate  toward  the  center.  This  geometry  (especially  in 
layered  sites)  then  does  not  lend  itself  particularly  well  to  fundamental 
study  of  the  shear  behavior  of  the  materials.  Comparisons  of  a  number  of 
calculations  of  explosive  events  with  experimental  data  indicate,  on  the 
other  hand,  that  the  shear  behavior  of  the  materials  plays  a  significant 
role  in  the  response  of  materials  to  explosive  loadings.  In  addition,  the 
importance  of  shear  waves  in  earthquake  problems  has  long  been  recognized 
and  has  precipitated  many  laboratory  studies  of  the  response  of  soils  to 
shear  loadings. 

These  laboratory  investigations  have  demonstrated,  conclusively,  that 
the  shear  modulus  of  soils  depends  strongly  on  the  shear  strain  to  which  the 
soil  is  subjected;  the  shear  modulus  decreases  significantly  as  the  maximum 
shear  strain  increases.  Therefore,  viable  insitu  testing  techniques  must  be 
capable  of  attaining  large  shear  strains.  These  laboratory  studies 
certainly  suggest  the  importance  of  the  shear  modulus  in  the  material 
response.  The  actual  determination  of  shear  properties  in  the  laboratory 
presents  significant  problems  associated  with  sample  disturbance  and  removal 
of  insitu  stresses  in  addition  to  those  associated  with  the  boundary 
conditions  imposed  by  the  laboratory  test  devices.  In  a  state-of-the-art 
address  to  the  ASCE  conference  on  "InSitu  Measurement  in  Soil  Properties" 
Wroth  (Ref.  11)  stated, 

"The  main  reason  for  the  need  for  insitu  measurement  of  the 
parameters  of  deformation  is  that  as  our  knowledge  of  the 
behavior  of  real  soil  expands,  so  does  our  appreciation  grow  for 
the  inadequacy  of  conventional  laboratory  testing.  The  marked 
consequence  of  the  inevitable  disturbance  that  is  caused  in  any 
soil  specimen,  however  carefully  it  is  sampled,  are  all  to 
evident." 
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Insitu  measurement  of  the  fundamental  properties  necessary  for  ground 
response  analysis,  therefore,  has  become  an  area  of  major  concern. 

Two  recent  papers  have  presented  an  excellent  review  of  the 

state-of-tne-art  in  the  in-situ  material  property  determination.  The  first 
of  these  by  Anderson  and  Espana  of  Fugro,  Inc.  (Ref.  12)  was  sponsored  by 
the  Electrical  Power  Research  Institute  and  focuses  on  the  entire  spectrum 
of  insitu  testing  techniques  from  the  earthquake  engineering  perspective. 
Therefore,  this  review  emphasized  the  measurement  of  shear  modulus.  The 

autnors  of  this  paper  divide  the  test  techniques  into  what  they  call  wave 

propagation  methods,  dynamic  system  response  methods,  and  direct/ load 
strength  metnods.  For  each  of  these  methods  the  basic  concept  is  discussed, 
the  methods  of  data  analysis  are  presented  and  the  operational  systems  are 
described.  A  number  of  proposed  test  techniques  are  also  presented  even 

though  they  have  not  yet  been  reduced  to  practice.  In  summarizing  their 
findings,  Anderson  and  Espana  stated, 

"A  comparison  of  the  three  general  classes  of  test  methods,  (wave 
propagation,  dynamic-system  response  and  direct/ load  strength) 
determined  that  no  single  existing  test  method  or  concept  for  a 
new  test  metnod  satisfies  all  the  requirements  for  a  high-strain 
amplitude,  dynamic  in-situ  property  determination." 

They  went  on  to  state,  in  conclusion,  that  although  a  number  of  techniques 
were  presently  available  for  making  dynamic  insitu  soil  property 
measurements  that  new  insitu  testing  methods  are  required  to  extend  the 
present  capability  and  to  varify  the  significance  of  various  earthquake 
loading  and  related  parameters  presently  thought  to  be  of  importance. 

The  second  report,  by  Lodde  (Ref.  13),  was  prepared  for  the  AFWL  and 
focused  on  only  those  techniques  for  evaluating  insitu,  high  amplitude, 
shear  behavior.  This  report  was  concerned  with  the  properties  required  for 
nuclear  weapons  effects  calculations.  It  concluded  that  since  the  shear 
stresses  which  could  be  propagated  were  limited  by  the  shear  strength  of  the 
material,  the  techniques  used  in  earthquake  studies  were  equally  applicable 
to  the  weapon  effect  problem.  Therefore,  after  a  review  of  the  existing 
laboratory  shear  tests,  this  report  reviewed  the  same  existing  insitu  test 
techniques  covered  by  Anderson  and  Espana.  Lodde  also  concluded  that  none 
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of  the  currently  existing  testing  procedures  are  completely  adequate  for  the 
problem  at  hand.  As  a  result,  the  current  state-of-the-art  still  involves  a 
combination  of  insitu  and  laboratory  testing.  In  all  proDability,  such  a 
condition  will  obtain  for  a  considerable  time  in  the  future.  However,  the 
corrections  currently  used  for  the  high  strain  portion  of  the  laDoratory 
data  to  account  for  insitu  properties  certainly  requires  evaluation. 

\ 

\ 

Both  of  these  studies  suggest  that  the  best  currently  available  single 
test  is  the  insitu  Impulse  Test  developed  by  Shannon  and  Wilson,  and 
Agbabian  Associates  for  the  Nuclear  Regulatory  Commission.  This  technique 
is  described  by  Miller,  et.  al.  in  Reference  14.  It  is  similar  in  nature  to 
the  cross-hole  seismic  test  but  involves  a  larger  magnitude  shear  input. 
This  input  is  generated  in  an  anchor  hole  approximately  10  inches  in 
diameter  by  impacting  an  anchor  seated  against  the  walls  of  the  bore  hole. 
This  impact  produces  a  vertically  polarized  shear  wave  which  propagates 
radially  and  is  detected  at  various  ranges  by  sensors  located  in  additional, 
smaller  diameter  bore  holes.  By  comparing  the  waveforms  as  a  function  of 
radius  from  the  center  of  the  loading,  the  characteri Stic  shear  wave  can  be 
identified  and  its  velocity  of  propagation  determined.  The  correspondi ng 
shear  strains  are  estimated  using  elastic  theory.  The  primary  disadvantage 
of  this  technique,  at  the  current  time,  is  tnat  the  largest  amplitude 
strains  are  generated  in  tne  immediate  vicinity  of  the  10  inch  diameter 
anchor  nole.  This  zone  has  been  significantly  disturbed  and  the  geostatic 
stresses  have  been  altered  by  the  drilling  of  the  anchor  hole.  Therefore, 
the  large  amplitude  strains  do  not  occur  in  an  undisturbed  material. 

The  major  disadvantage  of  current  techniques  which  attempt  to  generate 
high  amplitude  shear  strains  are  that  the  largest  strains  are  generated  in  a 
material  distrubeu  by  the  drilling  of  a  fairly  large  diameter  bore  hole  for 
placement  of  the  shear  generation  device.  The  CIST  also  would  suffer  from 
this  limitation  even  if  it  were  modified  to  generate  a  larger  shear  wave  by 
Dropagating  the  detonations  from  the  top  to  the  bottom  of  the  explosive 
hole. 


A  possible  solution  to  this  is  to  use  an  explosive  charge  configuration 
involving  three  charges  near  the  ground  surface.  The  detonation  sequence 
and  timing  of  the  charges  can  be  selected  so  that  the  explosive  energy  is 
directed  into  a  skewed  region,  thereby  creating  an  asymmetry  in  the 
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airblast.  Asymmetries  tend  to  generate  large  shear  waves.  A  schematic  plan 
and  profile  of  this  configuration  is  shown  in  Figure  5.1.  The  yields  of  the 
three  explosive  sources  could  be  similar  to  those  used  in  conventional 
refraction  surveys  but  could  be  varied  depending  upon  the  depth  of 
penetration  desired  and  the  nature  of  the  materials  involved.  The  relative 
location  of  the  charges  could  be  varied  such  that  the  pressure  when  the 
blast  wave  from  the  simultaneous  detonation  of  charges  1  and  la  impinges 
upon  charge  2,  it  woulo  be  approximately  equal  to  the  shear  resistance 
(strength)  of  the  soil.  Charge  2  would  be  detonated  just  as  the  blast  wave 
from  charges  1  and  la  reaches  it.  The  confinement  of  the  charge  2  effects 
by  the  airblast  from  charges  1  and  la  would  direct  the  energy  into  a  limited 
region.  This  will  create  an  asymmetry  which  will  generate  the  large 
amplitude  shear  waves  which  will  be  measured  by  instruments  located  in  bore 
holes  at  various  ranges  and  depths.  The  depths  and  ranges  of  measurement 
would  depend  upon  the  layering  of  the  site  and  the  levels  to  which  the  shear 
wave  measurements  were  desired. 

Since  no  large  diameter  bore  hole  is  required,  and  no  large  piece  of 
equipment  must  be  transported  from  site  to  site,  this  technique  should  be 
more  economical  than  any  of  the  existing  large  amplitude  techniques. 
Utilizing  this  explosive  technique,  shear  modulus  data  should  be  obtainable 
for  strains  up  to  those  corresponding  to  the  shear  failure  limits  of  the 
materials.  The  more  distant  measurements,  on  the  other  hand,  would  provide 
information  at  the  lower  strain  levels  so  that  a  comparison  with  existing 
low  strain  amplitude  in-situ  techniques  would  be  possible.  In  addition,  a 
much  larger  volume  of  material  would  be  tested  than  in  the  techniques 
utilizing  a  bore  hole  device  for  generation  of  the  shear  impulse.  For  this 
technique  the  cost  would  be  directly  proportional  to  the  number  of 
measurement  locations.  Therefore,  would  vary  with  the  type  of  application. 

The  snear  modulus  plays  an  important  role  in  the  response  of  soils  to 
dynamic,  complex  loadings.  Since  current  insitu  test  techniques  are  not 
capable  of  generating  the  large  shear  strains  required  for  constitutive 
modeling,  the  technique  discussed  provides  a  possible  significant 
enhancement  in  insitu  testing  capability. 
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Initial  Boundary  Loading  Conditions. 

Figure  5.1.  Explosive  Shear  Wave  Generator. 


5.2  Effective  Stress  Modeling 


5.2.1  Introduction 

Natural  soils  are  three  phase  media  consisting  of  solid  grains,  water 
and  air.  Most  constitutive  modeling  has  treated  soil  as  an  equivalent 
single  continuum  material.  It  is  natural  to  ask  if  improved  soil  material 
modeling  requires  a  more  detailed  and  specific  treatment  of  the  individual 
components  of  the  soil,  and  the  interaction  between  these  components.  The 
main  concern  is  over  the  effect  of  pore  fluid,  especially  the  water 
component,  on  the  behavior  of  the  soil.  There  are  three  main  areas  of 
concern: 

•  Effective  stress  and  soil  strength 

•  Pore  pressure  increase  and  potential  liquefaction 

•  Lubrication  effects  on  stiffness  and  strength  of  the  soil. 

The  following  paragraphs  discuss  some  of  the  main  points  associated 
with  effective  stress  and  liquefaction,  and  present  some  preliminary 
conclusions  with  regard  to  future  work  that  must  be  performed  in  order  to 
include  these  effects  in  material  models.  Lubrication  effects  are  not 
discussed  because  there  is  very  little  information  available  on  this 
subject.  It  is  believed,  however,  that  some  relatively  simple  experiments 
(dry  soil  versus  wetted  soil)  would  yield  insight  into  the  area.  Comments 
are  also  presented  on  the  computer  requirements,  in  terms  of  memory  and 
computing  time,  that  might  be  introduced  by  a  fuller  treatment  of  the 
individual  components  of  the  soil  media. 

5.2.2  Strength  and  Effective  Stress 

Of  tne  tnree  components  of  natural  soil,  only  the  soil  grain  matrix  has 
strength.  For  all  practical  purposes,  the  water  and  air  phases  can  only 
carry  hydrostatic  loads.  The  soil  grain  matrix  achieves  its  strength  from 
friction  a nd/or  intramolecular  forces  developed  between  the  particles.  Dry 
and  fully  drained  tests  on  natural  soils  show  a  dependence  of  strength  on 
tne  normal  stress  applied  to  the  sample.  When  a  soil  containing  pore  fluid 
is  loaded,  part  of  the  load  is  carried  in  the  grain  matrix  and  part  of  the 
load  is  carried  in  the  air-water  phase,  depending  on  the  relative 
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compressibility  of  the  individual  components  and  on  the  ability  of  the  pore 
fluid  to  flow.  If  fully  saturated,  the  majority  of  the  load  is  carried  in 
the  water  phase.  If  less  than  fully  saturated,  tnen  most  of  the  load  is 
carried  in  tne  grain  matrix.  Since  only  the  grain  matrix  has  strength  and 
since  tnat  strength  is  dependent  upon  the  normal  load  applied  to  the  grain 
matrix,  it  is  apparent  that  only  the  normal  stress  carried  through  the 
intergranular  contacts  influences  the  strength  of  the  natural  soil.  This  is 
known  as  the  effective  stress  principal  and  was  established  by  Terzagi  (Ref. 
15)  over  40  years  ago.  The  effective  stress  is  expressed  quite  simply  by 
the  equation 


a  =  c  -  u 


where  a  is  the  total  stress  acting  on  a  given  pla..e,  u  is  the  pore  pressure 
in  the  pore  fluid,  and  a  is  the  effective  stress  which  is  carried  in  the 
grain  matrix. 

Most  modeling  for  dynamic  purposes  (blast  and  earthquake)  has  not 
treated  tne  details  of  effective  stress.  This  is  because  it  has  normally 
been  assumea  that  the  dynamic  phenomena  of  interest  occurs  so  fast  that  the 
material  can  be  considered  to  oe  undrainea  during  the  dynamic  process.  As  a 
result,  flow  does  not  occur  and,  therefore,  the  effective  stress  in  the  soil 
does  not  cnange.  In  turn,  the  strength  does  not  change  during  the  dynamic 
process.  Hence,  total  stress  models  have  been  the  rule. 

In  recent  years,  however,  there  nas  been  increasing  concern  that 
effective  stress  should  be  treated  in  more  detail.  A  significant  model 
developed  witnin  the  blast  and  shock  community  is  that  of  Baladi  (Ref.  16 
and  Appendix  A. 1.2c).  Baladi  developed  an  effective  stress  model  for  ground 
shocx  computations  by  considering  two  models  for  the  material.  One  was  a 
total  stress  model  derived  from  undrained  tests  on  the  soil  material.  The 
second  was  an  effective  stress  model  derived  from  fully  drained  tests  on  the 
same  material.  It  was  assumed  that  lubrication  effects  on  the  stiffness  and 
strength  of  tne  material  was  not  important.  In  Baladi 's  model,  computations 
of  total  and  effective  stress  variables  are  made  using  the  total  stress 
moael  and  also  using  the  effective  stress  model.  Using  assumptions  of 
volumetric  strain  compatibility  and  effective  and  total  deviatoric  stress 
compatibility,  the  results  of  the  computations  from  the  two  models  are  used 


85 


APPLIED  PE/6RPCH  R/ZOCIRTE/.inC, 


AD-A111  901  APPLIED  RESEARCH  ASSOCIATES  INC  ALBUQUERQUE  NM  F/G  8/13 

FUNDAMENTAL  PROPERTIES  OF  SOILS  FOR  COMPLEX  DYNAMIC  LOADINGS, <U) 

SEP  81  M  C  OASS.  J  L  BRATTON*  C  J  HIGGINS  FR9620-80-C-0088 

UNCLASSIFIED  AF0SR-TR-82-0101  NL 

2  2 


to  determine  the  pore  pressure.  This  pore  pressure  as  a  function  of  time  is 
produced  as  part  of  the  output  of  the  ground  motion  calculation. 


A  detailed  review  of  Baladi's  model  indicates  that  it  really  does  not 
provide  much  more  additional  useful  information  beyond  that  available  from  a 
standard  total  stress  model.  This  is  because  the  ground  motions  computed  by 
his  method  would  be  exactly  the  same  as  those  computed  with  a  total  stress 
model  alone  (Baladi  has  not  made  such  a  comparison).  Also,  the  loads  or 
stresses  produced  on  a  structure  embedded  in  the  soil  would  be  exactly  the 
same  since  those  loads  are  total  stress  loads  rather  than  effective  stress 
loads. 


The  usefulness  of  Baladi's  model  is  limited  primarily  by  the  fact  that 
it  ooes  not  include  fluid  flow  effects,  which  are  potentially  important  even 
for  general  ground  motion  phenomena.  Baladi's  model  does  not  include  any 
flow  because  of  the  assumption  that  the  phenomena  of  interest  occurs  so  fast 
that  flow  is  not  important  during  the  time  of  interest.  As  stated  earlier, 
this  has  been  a  standard  assumption  throughout  the  blast  and  shock  community 
over  the  past  twenty  or  so  years.  However,  there  is  another  view  that  has 
never  Deen  studied  in  detail  regarding  the  importance  of  flow.  Ishihara 
(Ref.  17),  in -an  analysis  similar  to  that  of  Biot,  derived  equations  of 
motion  for  the  solid  and  fluid  phases  of  a  soil  by  representing  the  soil  as 
an  elastic  porous  matrix  saturated  with  water.  His  main  purpose  was  to 
derive  results  which  could  be  interpreted  in  terms  of  standard  soil 
mechanics  concepts  and  tests.  Two  kinds  of  compression  waves,  the  same  as 
those  found  by  Biot,  were  derived  in  a  revised  form.  Evaluation  of  the 
deformation  modulus  in  the  frequency  equation,  with  reference  to  specific 
soil  test  conditions,  indicated  that  the  first  kind  of  wave  travels  through 
the  solid-water  phase  system  causing  only  compression  of  the  solids  and 
fluid  but  without  any  change  in  the  pore  volume  due  to  flow.  The  second 
kind  of  wave  can  progress  only  when  pore  volume  change  takes  place  due  to 
flow,  and  is  associated  with  consolidation.  The  most  interesting  finding 
from  Ishihara' s  work  was  that  flow  appears  to  be  more  important  for  what  is 
normally  considered  to  be  high  frequency,  undrained  loadings. 

Ishihara  states  the  usual  assumptions  and  his  findings  very  clearly: 
"In  view  of  the  usual  hypothesis,  it  is  expected  that  the  wave  that  is 
encountered  at  higher  frequency  would  be  the  one  that  does  not  permit  the 
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drainage  to  take  place,  because  the  fluctuations  in  the  pore  pressure  due  to 
strain  occurs  so  rapidly  tnat  tnere  is  little  time  for  the  water  to  flow 
from  regions  with  higher  pore  pressure  to  those  of  lower  pore  pressure.  On 
the  contrary,  the  wave  at  lower  frequency  is  expected  to  produce  a  condition 
where  the  pore  water  is  allowed  to  move  through  the  soil  skeleton,  because 
tnere  may  be  sufficient  time  available  for  an  equilibrium  to  be  established 
for  the  fluctuation  of  the  pore  pressure.  However,  the  above  expectation  is 
in  complete  contrast  with  what  was  found  above.  Nevertheless,  the  situation 
can  also  be  considered  from  another  point  of  view.  It  can  be  anticipated 
that  at  higher  frequency  the  wavelength  is  short,  so  that,  although  there  is 
little  time  for  water  to  move,  the  distance  that  the  pore  water  must  travel 
is  short.  On  the  contrary,  at  lower  frequency,  although  there  is  enough 

time  for  tne  water  to  flow,  it  has  to  traverse  over  a  great  distance  because 
the  wavelength  is  now  large.  According  to  this  viewpoint,  the  stress 
condition  that  is  encountered  with  a  wave  of  lower  frequency  is  an  undrained 
condition,  and  the  stress  condition  that  is  realized  with  a  wave  of  higher 
frequency  is  a  drained  condition." 

Ishihara's  observations  from  his  theoretical  work  are  in  marked 
contrast  to  the  standard  assumption  based  upon  intuition  which  has  been 
employed  in  dealing  with  dynamic  problems  assuming  total  stress  models. 
Ishihara  analyzed  the  one  dimensional  consolidation  problem  under  a 
sinusoidally  varying  force  to  demonstrate  his  concept.  The  well  known 

consolidation  dimensionless  time  constant  Ty  (Ref.  15),  which  determines 
the  time  it  takes  to  reach  pore  pressure  equilibrium,  is  inversely 

proportional  to  the  square  of  the  distance  and  directly  proportional  to 

time.  Hence,  the  distance  seems  to  be  more  of  a  controlling  factor  for  pore 
pressure  equilibrium  or  water  flow  than  time,  supporting  Ishihara's 
hypothesis.  Although  Ishihara's  work  was  published  in  1967  it  has  not  been 
widely  evaluated  in  the  blast  and  shock  community.  Hence,  it  is  a  major 
conclusion  of  the  authors  here  that  his  work  be  given  a  thorough  evaluation 
in  next  years  effort  using  available  data  and  calculations. 

5.2.3  Grain  Matrix  Degradation  (Liquefaction) 

The  other  major  phenomena  which  is  dependent  upon  pore  fluid  behavior 
and  effective  stress  is  that  of  liquefaction.  Figure  5.2  shows  some  results 
from  a  computation  given  in  reference  16.  The  calculation  was  for  a 
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Figure  5.2  Comparison  of  total  and  effective  stress  paths  range 
15.6  meters.  (After  Ref.  16) 


radially  expanding  airblast  on  a  homogeneous  halfspace  of  saturated  sand. 
The  figure  shows  the  total  and  effective  stress  paths  at  three  different 
depths  at  a  range  of  15.6  m  from  the  origin.  One  of  the  most  interesting 
points  on  these  stress  paths  is  that  the  effective  stress  at  depths  of  1.6 
and  3.6  m  become  negative  during  unloading.  This  means,  in  a  physical 
sense,  that  the  grain  matrix  has  been  or  is  about  to  be  broken  down  or 
degraded.  Indeed,  this  is  the  onset  of  liquefaction  under  compressive 
loading  and  unloading.  Baladi  does  not  comment  on  this  behavior  at  all  in 
his  report.  However,  it  is  here  that  one  of  the  most  interesting  parts  of 
the  effective  stress  problem  begins  to  occur,  i.e.,  liquefaction. 

In  order  to  deal  with  liquefaction  and  post-liquefaction  stability  of 
soil,  a  model  must  be  able  to  handle  three  complex  problems: 

•  Air-water  compressibility 

•  Pore  pressure  changes 

•  Fluid  flow 

Air-water  compressibi 1 ity  is  a  major  consideration  because  many  soils, 
especially  in  regions  near  sea  coasts,  rivers,  waterways  or  where  water 
table  fluctuations  occur,  are  not  in  a  state  of  100  percent  saturation. 
Minute  amounts  of  air,  which  generally  are  impossible  to  measure  by  standard 
soil  mechanics  tests,  can  markedly  affect  the  compressibi ity  of  the 

air-water  pnase.  As  soon  as  the  compressibility  of  the  air-water  phase 
becomes  large  compared  with  the  compressibility  of  the  soil  grain  matrix, 
even  if  the  material  is  98  percent  saturated  and  above,  the  material  behaves 
essentially  as  a  soil  grain  matrix  only.  However,  depending  upon  the 

magnitude  of  the  load,  the  air  contained  in  the  air-water  phase  can  be 
driven  into  solution  and  the  material  can  become  100  percent  saturated. 

Tne  phenomena  by  which  these  changes  occur  is  extremely  complex 

involving  the  behavior  of  free  air  bubbles,  bubbles  adhered  to  soil  grains, 
surface  tension,  the  solubility  of  air  in  water  and  a  number  of  other 
factors  including  the  dynamic  behavior  of  these  parameters.  Air-water  phase 
compress i d i 7 ity  on  a  simplistic  scale  has  been  treated  by  Wood  (Ref.  18) 
and,  in  Russia,  by  Lyakhov  (Ref.  19).  A  recent  investigation  in  the 
Netherlands  by  Barends  (Ref.  20)  provides  a  much  more  sophisticated 

treatment  of  all  of  the  factors  involved  in  air-water  compressibility 
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including  most  of  the  factors  mentioned  above  but  not  on  a  dynamic  basis. 
Barends  work  indicates  tnat  the  air-water  phase  compressibi  i ity  does  not 
vary  continuously  but,  rather,  that  distinct  pressures  exist  beyond  which 
the  free  air  becomes  dissolved  quite  suddenly.  Bubble  collapse  causes  an 
unstable  increase  in  compressibility.  The  main  point  here  is  not  to  confuse 
the  issue  but  to  make  clear  that  air-water  compressibility  in  and  of  itself 
is  a  major  problem  in  the  development  of  multi-phase  models. 

To  illustrate  the  effect  of  small  amounts  of  water  upon  the 
compressibi 1 ity  of  the  air-water  phase,  the  pressure  volume  relationships 
for  air-water  mixtures  at  various  degrees  of  saturation  using  the  simpler 
model  of  the  Russians  (Ref.  19)  have  been  calculated.  The  results  are  shown 
in  Figure  5.3.  The  results  indicate  that  degrees  of  saturation  less  than 
99.99  percent  yield  a  relatively  large  range  of  pressures  over  which  the 
air-water  phase  shows  very  high  compressibility.  Indeed,  for  degrees  of 
saturation  below  about  95  percent,  the  material  can  be  considered 
essentially  void  of  water  except  for  any  lubrication  effects  that  might 
occur. 


The  second  major  concern  in  the  analysis  and  prediction  of 
liquefaction  is  the  development  of  a  model  which  will  predict  the  pore 
pressure  changes.  These  changes,  of  course,  are  due  to  changes  in  the  pore 

volume.  There  are  several  models  developed  within  the  earthquake  community 

which  deal  with  pore  pressure  changes  under  multiple  cycles  of  triaxial 
compression  or  simple  shear.  These  models  are  generally  restricted  to  these 
limited  stress  paths.  Because  of  the  simple  assumptions  in  the  models 

developed  in  the  earthquake  engineering  community,  they  cannot  be  used  for 
more  general  stress  paths,  particularly  those  related  to  blast  and  shock 
problems  where  compression  is  important.  On  the  other  hand,  there  are 
models  which  can  deal  with  liquefaction  under  the  influence  of  compressive 
loads.  These  models,  however,  cannot  predict  earthquake-type  behavior 

because  they  generally  lack  coupling  between  deviatoric  and  volumetric 
strain. 


The  earliest  work  that  has  come  to  our  attention  which  deals  with 
compression  (or  blast)-induced  liquefaction  is  that  performed  by  the 
Russians  (e.g.  Refs.  19  and  21).  This  work  was  done  in  the  1950' s  and 
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1960‘s.  More  recent  work  has  been  done  in  Europe  (e.g.  Ref.  22),  but  the 
European  models  are  very  much  like  tne  Russian  models. 

To  illustrate  what  might  occur  after  the  effective  stresses  go 
negative  in  Baladi's  problem  (fig.  5.2),  consider  the  following  simple 
illustration  based  upon  the  Russian  models.  Available  information  on 
explosion-induced  liquefaction  reveals  that  the  pore  pressure  behavior  is 
dependent  upon  the  same  parameters  which  are  important  for  earthquake 
phenomena,  i.e.,  air  content,  relative  density,  and  all  the  other  parameters 
which  influence  the  grain  skeleton  compressibility  (grain  size,  soil 
structure,  grain  roughness,  cementation,  etc.).  The  liquefaction  mechanism 
currently  accepted  by  the  Europeans  (e.g.  Ref.  22)  and  proposed  by  the 
Russians  in  the  1960’s  (Ref.  21)  is  based  upon  irrecoverable  strain  in  the 
soil  skeleton  upon  unloading  from  some  previous  volumetric  strain.  If  the 
fluid  cannot  flow  during  this  period,  residual  pore  pressures  develop. 

The  basic  mechanics  for  a  simple  model  for  composite  fluid-skeleton 
behavior  in  uniaxial  strain  is  given  in  Figure  5.4.  On  loading,  both  the 
fluid  and  the  skeleton  carry  stress.  The  total  stress  is  given  by 

a  *  a  (c)  +  p(^) 


where  o  *  total  stress 

T  =  stress  in  the  skeleton  (intergranular  stress) 
p  =  pore  pressure 

The  strain  in  the  soil  skeleton  is  the  total  strain  ,  while  the  strain  in 
the  pore  fluid  phase  is  given  by  e/n  where  n  is  the  porosity.  Because  the 
fluid  is  usually  much  stiffen  than  the  skeleton,  most  of  the  stress  on 
loading  is  carried  in  the  fluid.  Upon  unloading,  the  skeleton  exhibits 
irreversible  compaction  so  that  the  skeleton  stress  (effective  stress) 
reaches  zero  while  there  is  still  pore  pressure  in  the  fluid.  It  is  at  this 
point  that  the  soil  matrix  tends  to  break  down  or  liquefy. 

An  example  of  the  load-unload  behavior  for  a  soil  which  is  99  percent 
saturated  and  which  has  skeleton  compressibility  similar  to  that  of 
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Figure  5.4.  Concepts  for  a  simple  model  for  liquefaction  in  uniaxial  strain 


r 


McC  rmick  Ranch  sand  is  given  in  Figure  5.5.  The  composite  fluid-skeleton 
is  loaded  to  a  peak  total  stress  aDove  28  MPa.  Upon  unloading,  the 
skeleton,  because  of  residual  strain,  reaches  zero  intergranular  stress 
while  tne  fluid  phase  still  carries  over  13  MPa.  At  this  stage,  the  soil 
skeleton  breaks  down  and  further  unloading  occurs  in  a  heavy  fluid  composed 
of  grains  and  water.  In  a  fully  saturated  soil,  liquefaction  would  occur  at 
lower  stress  levels.  The  example  serves  to  illustrate  both  the  mechanism  of 
liquefaction  and  the  importance  of  the  air-water  phase  compressibility. 

It  is  the  relative  compressibi  1  ity  between  the  soil  skeleton  and  the 
air-water  phase  which  controls  whether  liquefaction  occurs  and  the  stress 
levels  needed  to  cause  it.  The  relative  compressibility  can  be  affected  by 
the  skeleton  stiffness  on  the  one  hand  or  the  air-water  phase 
compressibility  on  the  other.  Whenever  the  relative  compressibility  is  such 
that  the  skeleton  remains  elastic  (i.e.,  rebounds  fully  on  unloading)  or  is 
stiffer  than  the  air-water  phase,  liquefaction  will  not  occur. 

Tne  simple  model  for  liquefaction  in  a  load/unload  cycle  of  uniaxial 
strain  will  not  predict  liquefaction  under  shear  conditions,  just  as  simple 
shear  models  will  not  predict  liquefaction  in  the  uniaxial  case  just 
described.  In  order  to  provide  a  good  prediction  of  pore  pressure  buildup 
ana  tne  potential  for  liquefaction  under  general  stress  paths,  it  is 
necessary  to  have  a  model  for  the  grain  matrix  which  includes  a  good  model 
of  deviatoric  (shear)  strain-volumetric  strain  interaction  encompassing  the 
entire  range  of  potential  stress  paths  which  can  be  encountered  in  dynamic 
problems.  None  of  the  models  in  use  in  the  blast  and  shock  community  today 
include  deviatoric-volumetric  strain  interaction.  The  cap  model,  for 
example  supresses  volume  change  at  the  intersection  of  the  cap  and  the 
standard  yield  surface.  Models  commonly  need  for  earthquake  related 
problems  do  not  include  this  general  capability  either.  Critical  state 
models,  although  not  generally  in  use  for  dynamic  problems,  to  account  for 
volume  decrease  or  increase  which  is  known  to  physically  occur. 

5.2.4  Conclusions  and  Recommendations 

The  key  points  that  must  be  considered  for  effective  stress  modeling 
are,  first,  the  importance  of  flow  during  the  time  of  interest  and  its 
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Load  to  30.6  MPa 


Figure  5.5.  Illustration  of  the  mechanism  of  liquefaction  in  a  load 
unload  uniaxial  rain  cy' 
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influence  on  effective  stresses,  and  second,  the  ability  of  the  grain  matrix 
model  to  deal  with  deviatoric-volumetric  strain  coupling. 

If  flow  occurs  and  causes  a  reduction  in  pore  pressure  during  the  time 
of  interest,  then  those  properties  of  the  material  which  are  dependent  upon 
integranular  or  effective  stress,  especially  strength,  will  change.  In 
addition,  the  potential  for  liquefaction  will  also  be  affected.  As 
discussed  previously,  Ishihara's  work  suggests  that,  even  in  dynamic 
problens  associated  with  blast  and  shock,  flow  and  resulting  pore  pressure 
changes  may  well  be  important.  In  fact,  Ishihara  also  points  out  that  for 
earthquake-related  problems  where  frequencies  are  lower  flow  may  not  be 
important.  It  is  essential  that  this  position  be  evaluated. 

The  second  major  point  bears  on  the  ability  of  our  models  to  deal  with 
liquefaction  phenomena.  In  order  to  do  so,  it  is  necessary  to  have  a  good 
model  for  the  grain  matrix  which  includes  deviatoric-volumetric  strain 
coupling  over  the  entire  credible  set  of  stress  and  strain  paths  ranging 
from  simple  shear  through  hydrostatic  compression  in  multiple  loading 
cycles.  In  addition,  it  is  necessary  for  a  good  air-water  phase  model  to 
account  for  the  effect  of  air-water  phase  compressibility  on  liquefaction. 
Finally,  it  is  essential  to  have  flow  in  the  model  to  account  for  pore 
pressure  changes  which  might  affect  intergranular  stresses,  a.'  well  as  any 
pore  pressure  changes  that  would  occur  over  a  longer  period  of  time  after 
liquefaction  has  occured.  This  is  necessary  in  order  to  predict  ground 
motions  as  well  as  analyze  instability  of  large  bodies  of  soil  due  to 
liquefaction. 

Regarding  ground  motion  prediction.  Figure  5.2  from  Baladi  indicated 
that  the  intergranular  stress  on  unloading  went  to  zero  while  there  was 
still  pore  pressure  in  the  fluid.  At  this  point  it  is  expected,  physically, 
that  the  grain  matrix  will  break  down  and  the  resulting  mixture  will  be  a 
heavy  fluid  consisting  of  grains  and  water.  In  a  ground  motion  problem, 
whether  it  be  blast  and  shock  or  earthquake-related,  there  may  still  be 
important  waves  propagating  through  the  media.  These  waves  will  be 
propagating  through  a  heavy  fluid  and  not  through  a  solid,  as  we  normally 
assume.  Hence,  even  for  this  part  of  the  wave  propagation  problem,  it  is 
necessary  to  have  a  good  model  which  considers  flow  and  as  appropriate,  will 


96 


APPLIED  PE/ERPCH  A/jfOCIflTE/.mC. 


degenerate  to  a  heavy  fluid  under  proper  physical  conditions.  One  model 
tnac  we  are  aware  of  whicn  contains  flow  is  due  to  Hart  (Ref.  23). 

Hart's  model  couples  thermal,  mechanical  and  fluid  flow  behavior.  His 
treatment  of  flow  appears  to  be  reasonable.  The  model,  however,  does  not 
include  a  modern  representation  for  the  grain  matrix.  He  simply  uses  an 
elastic-perfectly  plastic  model.  In  order  to  make  it  applicable  to  most  of 
the  problems  of  interest  in  soils,  the  grain  matrix  model  must  be  updated  to 
include  nonlinearity  before  failure,  as  well  as  an  improved  model  of 
deviatoric-volumetric  strain  coupling.  In  addition.  Hart's  model  must  be 
updated  to  include  a  better  air-water  phase  component  model. 

In  summary,  there  is  evidence  to  suggest  that  it  is  important  to  model 
effective  stress  in  ground  motion  proDlems.  This  modeling  must  encompass 
more  than  just  the  modeling  of  effective  stress.  It  is  also  necessary  to 
model  flow.  The  model  must  recognize  the  transition  from  a  relatively  solid 
tri-component  medium  to  a  heavy  fluid  when  liquefaction  occurs,  subsequently 
track  waves  which  propagate  through  this  heavy  fluid,  and  finally,  track 
long  term  recovery  to  a  stable  state.  There  are  bits  and  pieces  available 
for  such  a  model  in  various  places.  All  of  the  bits  and  pieces  however, 
have  shortcomings.  A  program  to  overcome  these  shortcomings  will  be 
developed  in  the  next  phase  of  this  project. 

The  limited  experience  that  is  available  on  the  use  of  effective 
stress  models  in  computations  comes  from  Baladi's  work  and  from  some  work 
performed  at  ARA  using  Hart's  model.  This  work  indicates  that  there  are  not 
major  computer  requirements  associated  with  the  incorporation  of  effective 
stress  information  or  flow.  Core  requirements  might  increase  by  10  to  15 
percent  in  order  to  accommodate  additional  equation  of  state  variables,  as 
well  as  a  few  additional  response  variables.  Computation  time  might 
increase  on  the  order  of  20  to  25  percent.  It  is  believed  that  these 
additional  requirements  are  relatively  minor  when  compared  with  the  improved 
predictions  and  analysis  which  would  be  available  from  the  use  of  effective 
stress  models. 
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6.0  SUMMARY 


This  report  briefly  describes  the  status  of  our  research  into  the 
dynamic  uehavior  of  soils.  Witn  tne  initial  year's  effort  at  a  close,  we 
are  continuing  as  outlined  in  the  proposal  for  continued  study* 

The  soil  element  model  has  been  developed  and  will  serve  as  a  basic 
tool  for  studying  the  merits  and  applicability  of  various  constitutive 
laws.  This  computer  program  will  continue  to  improve  as  additional  features 
are  added. 

Several  two-dimensional  calculations  have  been  performed  in  order  to 
better  understand  insitu  soil  behavior.  The  effect  of  varying  material 
moaels  was  investigated  as  well  as  establishing  some  "best-fit"  models,  and 
this  effort  will  be  continued. 
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APPENDIX  B 


Definitions  and  Notation 


B.l  Coordinate  System 

For  the  soil  element  model  code,  a  consistant  set  of  coordinate  axes 
has  been  maintained  as  much  as  possible.  Z-axial,  R-radial,  and 
T-tangential,  have  been  used  as  subscripts. 


The  stress  and  strain  tensors  for  this  system  are  then: 


_°RR  °RZ  aRT“ 

“'RR 

GR7 

eRT“ 

[0]  - 

°RZ  °ZZ  °ZT 

DO  = 

eRZ 

ezz 

£ZT 

_  °RT  °ZT  aTT  _ 

-GRT 

eZT 

GTT. 

Summation  notation  is  occasionally  used  here  to  minimize  cumbersome 
equations. 

B.2  Invariants 

Stress  and  strain  invariants  are  useful  in  that  they  describe  the 
stress/strain  state  in  quantities  which  are  independent  of  coordinate  axes 
configuration.  Stress  invariants  Jg,  and  J3  are  defined  here, 

strain  invariants  Ij,  Ig *  and  I^  follow  similarly. 
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The  three  invariants  are  the  coefficients  of  the  cubic  equation: 


+  J2o  +  =  0 

where  o^,  o2,  and  o3  are  the  roots  of  the  cubic,  also  known  as 

principal  stresses,  and: 


J1  "  °RR  +  °ll  +  aTT 


ctrr  °rz 

aRR  aRT 

aZZ  aZT 

= 

+ 

+ 

aRZ  °ZZ 

aRT  °TT 

°ZT  °TT 

J3  “ 


aRR  °RZ  °RT 


°RZ  °ZZ  °ZT 


aRT  aZT  °TT 


B.3  Deviatoric  Stress/Strain 

The  stress  deviator  tensor  is  defined  as: 


(oRR  "  P) 

aRZ  aRT 

[si]  ■ 

°RZ 

( aZZ  "  P^  aZT 

°RT 

°ZT  (aTT  “  P) 

L— 

J 

where  P  =  pressure  * 

°RR  +  aZZ  +  aTT 

The  deviator  tensor  is  also  refered  to  as  [Si j]. 
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J2*anc*  J3;  are  the  invariants  of  the  stress  deviator  tensor.  A 
convenient  derivable  expression  in  terms  of  the  full  stress  tensor  is: 


Invariants  of  the  strain  deviator  tensor,  I^‘,  I2*  and  I 3 ‘ ,  follow 

similarly. 


B.4  Volume  Behavior 

As  stated  above  P  denotes  mean  confining  stress,  or  pressure. 
Volumetric  strain,  eKK,  is  defined  as: 


e  _  e  +  e  +e 
KK  RR  11  TT  " 


is  also  known  as  the  cubical  dilitation,  or  the  first  strain 

invariant. 

The  excess  compression,  y,  which  is  used  in  the  AFWl  model,  is  defined 
as: 


where  p  =  density,  and 
pQ  *  initial  density 
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B.5  Moduli i 


In  general,  only  three  elastic  modulii  are  used  for  defining  elastic 
oeiiavi or  in  this  study: 

K  =  Bulk  Modulus  =  P/ 

G  =  Shear  Modulus  =  (cZ2  -  aRR)/  (ez-  -  eRR) 

v  a  Poisson's  Ratio  =  ^  ~ 

2(3K  +  G) 


B.6  Sign  Convention 

The  SEM  uses  a  sign  convention  which  is  typical  of  finite 

difference/finite  element  codes  which  are  in  use  for  geotechnical 
engineering.  Tension  is  considered  positive  for  the  stress/strain  tensor 

quantities  and  compression  is  positive  when  describing  volume  behavior 
quantities,  such  as  mean  confining  stress.  This  is  somewhat  confusing  to 

follow,  but  necessary  if  the  mocels  which  are  developed  with  this  program 

are  to  be  later  implemented  in  larger  codes. 
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APPENDIX  C 


Sample  SEM  Exercises 


Six  examples  are  given  to  illustrate  use  of  the  Soil  Element  Model. 
Note  that  units  vary  (SI,  English,  etc.)  for  each  example.  Any  consistant 
set  of  units  may  be  used  in  the  Soil  Element  Model. 


Example  I 


Material:  McCormick  Ranch  Sand 

Model:  Elastic 

K  =  8375  psi 
G  =  3865  psi 

Loading:  a)  Initial  geostatic  pressure  »  0.0  psi 

Orained  isotropic  compression 
Load  to  20  psi 
Unload  to  0  psi 
Load  to  50  psi 

d)  Drained  Std.  Triaxial 

Load  to  1  percent  axial  strain  difference 
Unload  to  0  psi  stress  difference 
Load  to  3  percent  axial  strain  difference 
Unload  to  0  psi  stress  difference 
Load  to  7  percent  axial  strain  difference 


Example  II 

Material:  McCormick  Ranch  Sand 

Model:  Elastic  -  Plastic 

K  -  8375  psi 
G  .  3865  psi 
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Loading: 

Example  III 

Materi al : 

Mode  1 : 

Loadi ng: 


Example  IV 

Material: 

Mode  1 : 


C,  =  260  psi 
CQ  =  260  psi 
Cc  =  0.016  (psi)-1 
TCUT1  =  0.0  psi 
FC UT1  =  0.0  psi 

Same  as  Example  I 


McCormick  Ranch  Sand 


Pyke‘ s  1-D  Curve  Fit 
K  =  8375  psi 
•ty  =  210  psi 

Giax  -  50-000  »s1 


Initial  geostatic  pressure  -  0.0  psi 
a)  Drained  Isotropic  Pressure 
Load  to  50  psi 


b) 


Strain  Controlled  Pure 
strain  reversal) 

5  cycles  at  YRZ  = 
5  cycles  at  YRZ  " 
5  cycles  at  YR^  = 


Shear  (using 

1.0  percent 
0.4  percent 
2.0  percent 


McCormick  Ranch  Sand 

AFWL  Engineering 
pz  «  1.75  g/cc 
m  0.30 
Vy  m  0.30 
Ku  -  74,000  psi 


complete 
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=  5,400  psi 
K2  =  12,300  psi 
«2  =  74,000  psi 
=  74,000  psi 
Pj  =  0  psi 
P2  =  150  psi 
Pj  =  520  psi 
T  =  0  psi 
Y  =  0  psi 
S  =  0.54 
VM  =  250  psi 

Loading:  Initial  geostatic  pressure  =  0.0  psi 

a)  Drained  Isotropic  Compression 

Load  to  20  psi 
Unload  to  0  psi 
Load  to  50  psi 
Unloaa  to  0  psi 

b)  Uniaxial  Loading  (Drained) 

Load  to  Axial  Stress  =  100  psi 
Unload  to  Axial  Stress  »  0  psi 
Load  to  Axial  Stress  =  200  psi 
Unload  to  Axial  Stress  =  0  psi 
Load  to  Axial  Stress  =  250  psi 


Example  V 

Material:  McCormick  Ranch  Sand 

Model:  Cap  Model 

*  7.50  ksi 
K1S  =  0.94 

K2S  =  1.0  (ksi)'1 
Gj  *  40  ksi 

*  0.75 
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Loading: 


Example  VI 

Material : 

Mode  1 : 


G2  =  3.0  (ksi )-i 

a  =0.0 

C$  =  0.08  ksi 
8  =  0.214  (ksi)'1 

C  =  0.39  ksi 

R  =  2.5 

0  =  0.7  (ksi)-1 

W  =  0.066 


Initial  geostatic  pressure  =  0.0  ksi 
a)  Drained  Isotropic  Compression 
Load  to  0.010  ksi 
Unload  to  0.0  ksi 
Load  to  0.050  ksi 


b) 


Drained  Std.  Triaxial 

Load  to  1  percent  axial  strain  difference 

Unload  to  0  ksi  stress  difference 

Load  to  3  percent  axial  strain  difference 


Ocean  Mud  (Seabed  Deposits) 


Effective  Stress  Cap 

KiS  =  20.0  MPa 

Kls  =  0.999 

K2S  =  0.0002  (MPa)-1 

KiM  =  1875  MPa 

K1m  =  0.35 

K2M  =  0.05  (MPa)-1 

G1  .  0.60  MPa 

Gj  -  -1.0 

G2  .  200.0  (MPa)-1 
a  =  0.133336 
Cs  -  0.00866  MPa 
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8  =  0.0  (MPa)"1 

C  =  0.0  MPa 

K  =  4.0 

0  =  1.20  (MPa)"1 

W  =  0.42 

Loading:  Initial  geostatic  pressure  =  0.037  MPa 

a)  Drained  Isotropic  Compression 

Load  to  0.050  MPa 
Unload  to  0.037  MPa 

*  Load  to  0.086  MPa 

b)  Undrained  Std.  Triaxial  Test 

Load  to  2.0  percent  axial  strain 
difference 

Unload  to  0.0  MPa  stress  difference 
Load  to  4.0  percent  axial  strain 
difference 

*  Unload  to  0.0  MPa  stress  difference 
Load  to  8.0  percent  axial  strain 
difference 


The  following  figures  illustrate  material  response  for  the  six  examples 
as  outlined  above. 
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EXAMPLE  I  -  ELASTIC  MODEL 
MCCORMICK  RANCH  SfiND 


EPF  PRESS  VS.  VOLUMETRIC  STRAIN 


USOTAOFtC  C0HFRC5SI0N) 
o 


VOLUME  STRAIN  *10 


s 


us 
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EXAMPLE  I  ~  ELASTIC  MODEL 
MCCORMICK  R  FINCH  SfiKD 

STRESS  RATIO  VS. AXIAL  STRAIN 


ISTRICWD  TCIflxlRL  TEST) 


VOLUME  STRAIN 
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EXflMFLE  I  -  ELASTIC  MODEL 
MCCORMICK  RANCH  SAND 
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EXAMPLE  II  -  ELASTIC-PLASTIC  MODEL 
MCCORMICK  RANCH  SAND 

EFF  PRESS  VS.  VOLUMETRIC  STRAIN 

I  ISOTROPIC  COMPRESSION) 
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EXAMPLE  II  -  ELASTIC-PLRSTIC  MODEL 
MCCORMICK  RANCH  5 FIND 
STRESS  DIFFERENCE  VS.  EFF  PRESS 


O  ISIflMJRRD  TRIRXIftL  TCSTJ 


EFF  PRESSURE  (PSD 
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II  -  ELRSTIC-PLfiSTIC  MODEL 
MCCORMICK  RANCH  SAND 
DIPF  VS.  RXIflL  STRAIN 


ISTRMDffiD  TR[ flXIRL  TEST) 
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EXAMPLE  II  -  ELASTIC-PLASTIC  MODEL 
MCCORMICK  RANCH  SAND 

CFF  PRESS  V5.  VOLUMETRIC  STRRIN 

ISTfiHOfRO  IRJftXIM.  TEST) 
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VOLUME  STRAIN 


131 


EXAMPLE  II  -  ELASTIC- PLASTIC  MODEL 
MCCORMICK  RANCH  SAND 

AXIAL  STRAIN  V5  VOLUME  STRAIN 

ISTAMDfRD  TRlflXIRL  TCST) 


AXIAL  STRAIN 
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EXAMPLE  II  -  ELASTIC-PLRSTIC  MODEL 
MCCORMICK  RfiNCH  SRND 

STRESS  RATIO  VS.RXIAL  STRAIN 

ISTRHDflRD  TRIflXIRL  TEST) 


AXIAL  STRAIN 
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EXAMPLE  II  -  ELA5TIC-PLRSTIC  MODEL 
MCCORMICK  RANCH  SfiND 

EPF  RXL  STRESS  VS.  RXIRLSTRR3N 

ISTfiMOFfiD  TRl fiXIRL  TEST) 


RXJRL  STRRIN 
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EXAMPLE  III  -  1-D  CURVE  FIT  (PTKE) 
MCCORMICK  RANCH  SAND 

EFF  PRESS  VS.  VOLUMETRIC  strain 

1ISDTR0PIC  COHPRESSIDM) 
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EXAMPLE  III  ~  1-D  CURVE  FIT  (PYKE) 
MCCORMICK  RRNCH  SRND 

SHEAR  STRESS  VS  .SHEAR  STRAIN 

tPUFE  SHERR  TESY) 
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EXAMPLE:  IV  -  FITHL  MODEL 
MCCORMICK  RANCH  5RN0 

EFF  PRESS  VS.  VOLUMETRIC  STRRIN 


I  ISOTROPIC  CONPRCSSIDN) 
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example:  iv  -  Arm  model 
MCCORMICK  RANCH  SAND 
EFF  PRESS  VS.  VOLUMETRIC  STRRIN 


(UNIflXIFL  7 ESI) 


EXAMPLE:  IV  -  APWL  MODEL 
ICODRMICK  RANCH  SAND 

STRESS  DIFP  VS.  RXIRL  STRRIN 


EXAMPLE  iv  -  AFWL  MODEL 
MCCORMICK  RANCH  SAND 

STRESS  OIFF  VS.  STRAIN  DIFF 


EXAMPLE  IV  -  RTML  MODEL 
MCCORMICK  RRNCH  SRND 

VOLUMETRIC  STRRIN  VS  STRRIN  DIFF 

tUNIRXIfL  IEST) 


VOLUME  STRAIN 
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EXAMPLE  IV  -  RTWL  MODEL 
MCCORMICK  RANCH  SAND 

AXIAL  STRRIN  VS  VOLUME  STRRIN 

(UNIRXIPL  1ESI) 
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EXAMPLE  IV  -  FirNL  MODEL 
MCCORMICK  RftNCH  SftND 

STRESS  RATIO  VS. AXIAL  STRAIN 


EXRMPLl  >V  ~  RTML  MODEL 
MCCORMICK  RANCH  SAND 

EFF  flXL  STRESS  VS.  RXlflLSTRRIN 


EXAMPLE  V  -  CAP  MODEL 

MCCORMICK  RANCH  SAND 

STRESS  DIFFERENCE  VS.  EFF  PRESS 

ISTBMOFRD  TRlflXIRL  TCST) 


MCCORMICK  RANCH  SAND 

STRESS  DIPP  VS.  RXIRL  STRRIN 

ISTflNOFRD  TRIRXIH! .  TCST) 


EXAMPLE  V  -  CAP  MODEL 

MCCORMICK  RANCH  SAND 

STRESS  DIFF  VS.  STRAIN  DIFF 

ISTRMDFRO  TRlflXIRL  TCST) 


EXAMPLE  V  -  rp-  '  .31  L_ 

MDCDRMI RANCH  SAND 

Err  PRESS  vs.  volumetric  strain 
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EXAMPLE  V 


-  CP.r  .3D-:. 

MDCDRMI  i;\  RANCH  SAND 

VOLUMETRIC  STRAIN  V5  STRFUN  QIFf 

ISTRHDPRO  TRlflKIRL  TEST) 
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flXIRL  STRRIN 


EXAMPLE  V  -  CAP  MODEL 

MCCORMICK  RANCH  SAND 

STRESS  RATIO  VS.RXIRL  STRAIN 

ISTRNDfftD  TRlflXIftL  TESTi 


EXAMPLE  V  -  CAP  MODEL 

MCCORMICK  RANCH  SAND 

EFT  aXL  STRESS  VS.  AXIflLSTRfllN 

I  STANDARD  TRfflXIftL  TESH 
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EXAMPLE  V  -  CAP  MODEL 

MCCORMICK'  RANCH  SAND 
EFP  RXL  STRESS  VS.  RXIRLSTRRIN 

^  I3THHCFRD  TRIflXIRL  TEST* 
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AXIAL  STRAIN 


153 


EXAMPLE  VI  ~  EFFECTIVE  STRESS  CfiP 
OCEAN  HUD 

EPP  PRESS  VS.  VOLUMETREC  STRAIN 


I ISDTRCPIC  COMPRESSION) 
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EXAMPLE  VI  -  EFFECTIVE  STRESS  CRP 
OCEAN  MUD 

TOTAL  RNO  EPF  STRESS  PATHS 

ISTfiMDFRD  TRlflXIFlL  TEST) 
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EXAMPLE  VI  -  EFFECTIVE  STRESS  CRP 
OCEAN  MUD 

STRESS  DIFF  VS.  RXIRL  STRAIN 

ISTAMDnftD  1R I AXIAL  TEST) 


RXIRL  STRREN 


157 


EXAMPLE  VI  -  EFFECTIVE  STRESS  CRF 
OCERN  MUD 

STRESS  DIFF  VS.  STRRIN  DIFF 

ISTRHDFRD  TRIflXIflL  TEST) 
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EXftfIPLE  VI  -  EFFECTIVE  STRESS  CAP 
OCEAN  MUD 

AXIfiL  STRRIN  VS  VOLUME  STRAIN 
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EXAMPLE  VI  ~  EFFECTIVE  STRESS  CRP 
OCEAN  HUD 

PORE  PRESSURE  VS.  AXIAL  STRAIN 


ISTBICFRD  TRtfiXIBL  TEST) 
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EXAMPLE  VI  -  EFFECTIVE  STRESS  CAP 
OCEAN  MUD 


STRESS  RATIO  VS. AXIAL  STRAIN 

ISTRHQFRD  TRtflXIRL  TEST) 
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EXAMPLE  ' 

EPP  nXL  s 


-  EFFECTIVE  STRESS  CAP 
''•DEPlN  MUD 

.SS  VS.  FIXIflLSTRRIN 


ISTRHDRRD  TRtRXIRL  TEST) 


11.13.27  KM  2  MOT,',  1981  JG8»HfVSfMJ  ,  (irWL  CRf)  DUSFtfl '£R  9.D 


EXFIMFLE  VI  -  EFFECTIVE  STRESS  CAP 
OCEAN  HUD 

TOT  RXL  STRESS  VS.  RXIflLSTRfllN 

ISIRMDFRD  IRIfiXIRL  TCST) 


RXIRL  STRRIN 


APPENDIX  D 


Professional  Personnel  and  Interactions 


D.l  Professional  Personnel 


The  following  people  were  associated  with  this  research  effort: 


1.  Jimmie  L.  Bratton 

Principal,  Applied  Research  Associates,  Inc. 

2.  William  C.  Dass 

Research  Engineer,  Applied  Research  Associates,  Inc. 

3.  Dr.  Cornelius  J.  Higgins 

Principal,  Applied  Research  Associates,  Inc. 

4.  Peter  Dzwi lewski 

Senior  Enyineer,  Science  Applications,  Inc. 

(presently.  Senior  Engineer,  Applied  Research  Associates,  Inc.) 

5.  David  K.  Rudeen 

Scientist,  Science  Applications,  Inc. 


6.  Dr.  Yoshi  Moriwaki 

Senior  Engineer,  Science  Applications,  Inc. 
(presently  with  Woodward-Clyde  Associates) 


D.2  Interactions 


A  symposium  entitled  "Implementations  of  Computer  Proceedures  and 
Stress-Strain  Laws  in  Geotechnical  Engineering"  was  held  in  Chicago, 
Illinois  on  August  3-6,  1981,  cosponsored  by  the  Illinois  Institute  of 
Technology  and  the  Virginia  Polytechnical  Institute.  The  conference  was 
attended  by  Mssrs.  Dass  and  Dzwi lewski  of  ARA,  and  provided  additional 
insight  into  the  present  status  of  research  on  constitutive  laws. 

In  addition,  ARA  has  submitted  two  abstracts  for  (see  below)  acceptance 
to  the  "International  Conference  on  Constitutive  Laws  for  Engineering 
Materials:  Theory  and  Application",  to  be  held  in  January  1983  at  the 
University  of  Arizona.  ARA  will  continue  to  prepare  papers  for  the 
appropriate  conferences  and/or  publications. 
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COMPARISON  OF  CONSTITUTIVE  MODELS 
SUBJECTED  TO  VARIED  STRESS-STRAIN  PATHS 

Willi  am  C.  Dass* 


ABSTRACT 


Many  levels  of  constitutive  models  are  available  today  to  the 
practitioner,  from  very  complex  to  simple  linear-elastic.  As  a  result, 
there  is  often  confusion  when  the  time  comes  to  choose  a  model  which  will 
give  the  right  combination  of  simplicity  (and  therefore  cost-effectiveness) 
and  accuracy. 

This  paper  describes  a  small  computer  program  which  has  been  developed 
as  a  tool  for  studying  material  constitutive  models.  The  program  allows 
comparisons  of  model  behavior,  and  with  it  one  can  perform  parametric 
studies  to  determine  the  influence  of  modeT  components.  It  is  intended  that 
the  "Soil  Element  Model",  as  it  is  called,  be  an  aid  in  model  development 
and  also  in  choosing  a  model  wnich  best  suits  a  particular  boundary  value 
problem. 

A  study  is  presented  which  illustrates  the  use  of  this  code  to  compare 
the  ability  of  several  material  models  to  replicate  laboratory  test  data. 
The  test  oata  (both  static  and  dynamic)  was  taken  from  previous  efforts  and 
some  data  was  generated  specifically  for  this  purpose.  The  study 
concentrates  on  parametric  effects  and  on  isolating  areas  for  model 
improvement  and  development.  In  addition,  complex  loading  paths,  typical  of 
those  which  are  induced  by  blast  loadings,  are  imposed  on  the  models  and 
some  comparisons  are  drawn. 
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THE  ROLE  OF  IN-SITU  TESTING  IN  DEVELOPMENT 
OF  CONSTITUTIVE  LAWS  FOR  ENGINEERING  MATERIALS 


Jimmie  L.  Bratton 


1 


ABSTRACT 


Theoretical  development  of  Constitutive  Laws  for  engineering  materials 
begins  with  a  set  of  characteristics  assumed  to  be  representative  of  the 
benavior  of  the  material  under  boundary  conditions  approximating  those  of 
the  engineering  problem  of  i rarest.  Generally  these  characteristics  are 
developed  from  laboratory  tests  of  the  material.  Laboratory  tests  of 
geologic  materials,  however,  are  subject  to  many  limitations.  The  sampling 
and  specimen  preparation  procedures  can  lead  to  both  quantitative  and 
qualitative  changes  in  the  material  behavior.  In  addition,  the  boundary 
conditions  which  can  be  applied  and  measurements  which  can  be  made  are 
limited  and  not  always  the  most  indicative  of  the  material  behavior  for  the 
conditions  of  interest  or  sufficient  to  allow  the  generalization  of  the 
measured  response  necessary  for  development  of  a  constitutive  model. 

In-situ  testing  has  been  shown  to  provide  a  necessary  and  important 
adjunct  to  laboratory  testing  for  soils  and  rock.  This  is  particularly  true 
when  the  loadings  of  interest  are  complex  and  dynamic.  In-situ  testing 
techniques  are  Driefly  reviewed  and  their  application  to  the  establishment 
of  the  material  response  characteristics  for  constitutive  modeling  is 
discussed.  The  use  of  these  tests  for  evaluating  models  is  also  reviewed 
and  specific  examples  are  presented.  Characteristics  of  the  material 
response  which  were  not  identified  in  laboratory  testing  are  presented  for 
specific  cases  and  the  role  of  in-situ  testing  in  the  development  of 
constitutive  laws  for  engineering  materials  is  illustrated. 
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